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Abstract 

This  thesis  was  concerned  with  the  development  of  a  computer  model 
of  a  passive  synthetic  aperture  imaging  system.  The  research  was 
divided  into  three  parts.  They  were  (1)  applying  an  understanding  of 
partial  coherence  theory  and  its  relationship  to  the  impulse  response  of 
the  system,  (2)  developing  the  computer  model,  and  (3)  exercising  the 
computer  model  to  perform  a  sensitivity  analysis. 

The  system  modeled  consisted  of  two  lenses  mounted  on  a  movable 
platform.  The  lenses  were  separated  by  a  fixed  distance  and  travelled 
in  a  direction  parallel  to  this  separation.  The  coherence  of  radiation 
present  at  each  lens  emanating  from  a  real  source  was  measured  yielding 
the  Fourier  transform  of  the  source  intensity  distribution  according  to 
the  van  Cittert-Zernike  theorem  (2:510).  The  transform  was  then  multi¬ 
plied  by  an  effective  aperture  (obtained  from  the  motion  and  position  of 
the  lenses  relative  to  the  source) .  An  inverse  Fourier  transform  was 
then  applied  to  this  result  yielding  the  image.  This  is  the  process 
modeled  by  the  computer. 

The  results  Indicated  that  new  means  of  image  Interpretation 
must  be  developed  in  order  to  make  the  results  useful.  This  is  due  to 
the  fact  that  the  system  behaves  much  like  a  high  pass  filter  and  the 
image  is  edge  enhanced  and  not  a  scaled  version  of  the  geometric  image. 


I.  Introduction 


The  goal  of  this  research  was  to  perform  a  sensitivity  analysis  of 
the  imaging  performance  of  a  passive  interferometric  imaging  system.  A 
hypothetical  system  consisting  of  two  lenses  that  are  physically 
connected  yet  separated  by  a  fixed  distance  was  determined  to  be  the 
simplest  case.  Such  a  system  could  be  mounted  on  a  movable  platform. 

The  system  samples  the  partially  coherent  radiation  emanating  from  a 
source  as  its  f ield-of-view  travels  across  the  source. 

A  cross-correlation  is  performed  between  the  radiation  fields 
present  at  each  lens  at  predetermined  intervals  of  time.  This  results 
in  a  set  of  discrete  samples  of  the  Fourier  transform  of  the  source 
radiant  intensity  distribution  being  measured  (4:2-1).  The  source 
radiant  intensity  distribution  may  then  be  found  by  taking  the  inverse 
Fourier  transform.  This  is  a  direct  application  of  the  van 
Cittert-Zernike  theorem  (2:510). 

The  technique  of  passive  interferometric  imaging  is  not  new. 

Radio  astronomers  have  used  it  for  quite  some  time  to  obtain  the  angular 
diameter  and  brightness  distributions  of  celestial  bodies  (13:2115). 
Efforts  to  Implement  this  technique  at  optical  or  infrared  frequencies 
have  yielded  promising  results  (14:1).  Current  synthetic  aperture 
systems  require  the  transmission  of  a  coherent  signal  to  obtain  an  image 
(4:2-1).  Attention  is  now  being  given  to  passive  interferometric 
techniques  in  the  infrared  region  because  this  technique  does  not 
require  the  use  of  an  active  illumination  which  could  reveal  the 
detector's  presence.  High  resolution  has  been  obtained  in  the  radio 
frequency  region  (13:2114)  and  it  is  desired  to  see  if  this  can  be 
duplicated  at  infrared  frequencies. 
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Goal  of  the  Thesis 


The  goal  of  this  thesis  was  to  develop  the  necessary  background  and 
methods  to  perform  a  sensitivity  ahalysls  and  to  use  this  information  to 
form  a  computer  model  of  the  system.  This  was  done  by  constructing  a 
model  of  the  overall  system.  The  computer  model  was  then  exercised  to 
see  how  varying  operating  conditions  and  parameters  that  the  imaging 
performance  is  dependent  upon  affected  system  performance.  These 
results  were  compared  to  the  results  obtained  under  what  were  defined  as 
ideal  operating  conditions  to  see  what  the  effects  were. 

The  thesis  was  divided  into  three  phases  in  order  to  meet  these 
broad  objectives.  They  were  (1)  applying  an  understanding  of  partial 
coherence  theory  and  its  relationship  to  the  impulse  response  of  the 
system,  (2)  developing  the  computer  model,  and  (3)  exercising  the 
computer  model  to  perform  the  sensitivity  analysis.  The  details  of 
these  phases  are  enumerated  below. 

Phase  One.  Understanding  partial  coherence  theory  revolves  around 
an  understanding  of  the  propagation  of  the  mutual  intensity  function 
(1:31).  This  is  the  quantity  measured  by  the  system.  A  normalized 
version  of  the  mutual  intensity  function  is  related  to  the  Fourier 
transform  of  the  source  intensity  distribution.  The  image  is  then  found 
by  taking  the  inverse  Fourier  transform  of  the  detected  normalized 
mutual  intensity  function.  Understanding  partial  coherence  and  the 
mutual  intensity  function  and  how  it  is  measured  are  therefore  the  first 
important  steps  in  analyzing  the  problem. 


Understanding  the  contribution  of  the  van  Cittert-Zernike  theorem 

is  the  next  step  in  obtaining  the  impulse  response  of  the  system.  As  is 

stated  in  this  theorem,  it  is  assumed  that  the  source  is  spatially 

Incoherent  and  emits  quasi-monochromatic  light.  The  system  behaves 

somewhat  like  a  coherent  imaging  system  despite  the  source  being 

incoherent.  This  is  illustrated  in  the  following  paragraph. 

It  can  be  shown  (7:110-113)  that  for  a  given  object  amplitude 

distribution  u(X  ),  the  resultant  coherent  image  amplitude  distribution 
o 

u(X^)  is  given  by 

u(X1)  -  u(Xq)  *  h(X)  (1.1) 

where  *  denotes  a  convolution  and  h(X)  is  the  amplitude  impulse  response 
of  the  imaging  system. 

The  amplitude  impulse  response  h(X)  is  the  Fourier  transform  of  the 
pupil  function  p(Adx)  (7:105)  where  A  is  the  wavelength  and  d  is  the 
distance  from  the  exit  pupil  of  the  optical  system  to  the  image  plane. 
Taking  the  Fourier  transform  of  both  sides  of  Eq  (1.1)  yields  the  linear 
system  equation 

U^)  =*  U(fQ)P(Adf)  (1.2) 

where  U(f.),  U(f  ),  and  P(Adf)  are  the  Fourier  transforms  of  u(X.), 
l  o  l 

u(Xq),  and  p(Adx)  respectively. 

This  shows  that  the  Fourier  transform  of  the  image  amplitude 
distribution  is  directly  proportional  to  the  Fourier  transform  of  the 
object  amplitude  distribution.  The  proportionality  constant  is  the 
scaled  pupil  function.  As  noted  earlier,  the  van  Cittert-Zernike 
theorem  states  that  the  mutual  intensity  function  is  the  Fourier 


transform  of  the  source  intensity  distribution.  Therefore,  applying  the 
van  Cittert-Zernike  theorem  to  a  given  source  intensity  distribution 
yields  the  input  to  the  linear  system  denoted  in  Gq  (1.2).  The  transfer 
function,  P(Adf),  is  then  the  only  parameter  needed  to  determine  the 
image  amplitude  distribution.  Once  P(Adf)  has  been  determined,  U(f^) 
can  be  calculated  for  any  source  amplitude  distribution  (via  Eq  (1.2)) 
and  u(X^)  can  be  found  by  taking  the  inverse  Fourier  transform  of  U(f^). 

This  first  phase  identified  the  parameters  and  conditions  which  may 
effect  the  impulse  response  and  system  performance.  The  parameters 
initially  identified  were  the  distance  between  the  lenses,  aperture 
size,  aperture  fill,  and  frequency  spacing  or  sample  spacing. 

Phase  Two.  The  second  phase  was  the  development  and  testing  of  the 
computer  model.  The  scenario  upon  which  the  model  was  based  is  as 
follows.  A  collection  platform  has  two  lenses  mounted  on  it.  These 
lenses  are  identical  and  are  separated  by  a  fixed  distance.  The  path 
travelled  by  the  collector  is  parallel  to  this  separation.  The  slant 
range  to  the  target  is  large  enough  such  that  it  lies  in  the  same  plane 
as  the  target  to  be  Imaged.  The  collector  traverses  an  angle  A0  as  it 
moves  past  the  target.  This  scenario  is  depicted  in  Figure  1.1. 

The  collector  ideally  would  be  able  to  gather  data  along  A6  equal 
from  0  to  180  degrees.  The  mutual  intensity  function  would  be  sampled 
all  along  this  interval  with  a  diffraction  grating  providing  the  means 
of  collecting  several  frequencies  each  time  a  sample  is  taken.  A  point 
source  was  considered  the  simplest  case  since  this  would  provide  the 


Phase  Three.  The  last  phase  consisted  of  exercising  the  computer 
model  to  obtain  point  spread  functions  of  the  system  as  well  as  images 
of  certain  simple  objects.  This  was  done  by  altering  the  ideal  model  of 
the  second  phase  with  changes  to  parameters  identified  in  the  first 
phase.  The  goal  was  to  determine  the  system's  imaging  performance  under 
realistic  operating  conditions. 


II.  THEORY  OF  OPERATION 


This  chapter  will  better  describe  and  more  fully  explain  the 
theoretical  basis  of  the  optical  system  introduced  in  the  previous 
chapter.  It  is  assumed  that  the  reader  has  an  understanding  of 
geometrical  and  Fourier  optics.  The  chapter  is  divided  into  three 
sections.  These  are  a  physical  and  conceptual  description  of  the 
optical  system  being  modeled,  the  theoretical  basis  of  operation,  and 
linking  the  theory  and  the  idea  of  a  coherent  imaging  system  to  the 
optical  system  being  modeled. 

The  Optical  System 

The  optical  system  consists  of  two  main  components.  One  is  a  pair 
of  lenses  separated  by  a  distance  d.  The  other  is  a  linear  array  of 
detectors.  The  goal  is  to  produce  an  image  of  a  thermal  source  which 
emits  a  randomly  fluctuating  field  by  measuring  the  complex  degree  of 
partial  coherence  of  the  radiation  field  present  at  the  lenses.  The 
system  can  be  conceived  as  one  mounted  on  a  movable  platform  that 
enables  the  system  to  rotate  around  the  scene  of  interest.  This  is 
conceptually  the  same  as  letting  the  system  remain  stationary  and 
allowing  the  scene  to  rotate  below  as  shown  in  Fig.  2.1.  Appendix  A 
relates  this  simple  geometry  to  the  more  complicated  geometry  of  the 
moving  lenses.  Enough  data  must  be  collected  so  that  an  image  of 
sufficient  resolution  and  quality  may  be  obtained  by  taking  the 
two-dimensional  inverse  Fourier  transform  of  the  collected  data. 


Fig  2.1.  System  Geometry 

The  scene  A  Is  assumed  to  be  within  the  field  of  view  of  the 


system,  and  L 2  are  the  two  lenses  of  the  system,  d  is  the  distance  by 
which  the  lenses  are  separated,  and  R2  the  distances  from  the  center 
of  the  field  of  view  to  and  respectively. 

Theoretical  Basis  of  Operation 

The  simplest  point  from  which  to  begin  is  Young's  two  slit 
experiment  (1:7-11)  which  will  be  used  to  introduce  the  mutual  intensity 
function  and  the  complex  degree  of  coherence.  Consider  the 
one-dimensional  setup  illustrated  in  Figure  2.2,  where  S  is  an  extended 
polychromatic  source  in  the  X  plane. 


►/ 

t; 


Fig  2.2.  Young' 8  Two  Slit  Experiment 
S^  and  S^  are  identical  slits  separated  by  a  distance  d  in  the  A 
plane  which  is  otherwise  opaque  and  is  parallel  to  and  a  distance  D  from 
the  X  plane.  Q  is  a  point  in  the  B  plane  which  is  parallel  to  and  a 
distance  Z  from  the  A  plane.  and  are  the  distances  from  and 
respectively  to  Q.  A  scalar  treatment  is  sufficient  here  because  the 
angles  involved  are  assumed  to  satisfy  the  small  angle  approximation 
(sin  0  s  6). 

The  source  S  emits  a  complex  incoherent  light  disturbance  E(P,t), 
This  disturbance  propagates  from  the  source  to  the  plane  A  according  to 
the  wave  equation  (1:9) 

(Y2)E  -  (1/c)2  d^E  (2.1) 

dt 


where  c  is  the  speed  of  light.  The  amplitudes  of  the  disturbance  at  the 


openings  S1  and  S2  are  then  denoted  E^(t)  and  E2(t)  respectively.  The 
total  disturbance  at  Q  is  then 

EQ(t)  -  B1(t-R1/c)  +  E2(t-R2/c)  (2.2) 

where  R^/c  and  R2/c  represent  the  time  delays  of  E^(t)  and  E2(t) 
respectively  in  propagating  to  Q.  Therefore,  let  R^/c  *  t^  and 
R2/c  ■  t2*  Rewriting  Eq  (2.2)  with  these  changes  yields 

EQ(t)  -  E1(t-t1)  +  E2(t-t2)  (2.3) 

The  irradiance  at  Q,  IQ,  is  necessarily  a  long  time  average  of  EQ(t). 
This  is  due  to  the  fact  that  the  frequency  of  the  radiation  field  being 
sampled  (EQ(t))  far  exceeds  the  capability  of  the  detectors  employed  to 
detect  each  individual  oscillation.  IQ  is  defined  as 

IQ  -  <(E1(t-t1)  ♦  E2(t-t2)  (E^t-tp  +  E2(t-t2))*>  (2. A) 

where  <.,.>  indicates  the  long  time  average  of  the  quantity  they 
enclose;  i.e. 

T 

<I(t)>  -  Lim  (1/T)  /I(t)dt  (2.5) 

-0 


Mutual  Coherence  Function.  Carrying  out  the  operations  indicated  in  Eq 


$ 


r-4»r*  t*.  «r-  ^  .  *i.  ^  .  - .  •-■-^  »  -  -  -.- 


( t 2 )  between  E,  and  by  t  (since  E^  and  E^  are  assumed  to  be 
stationary  fields),  Eq  (2.6)  can  be  rewritten  as 

IQ  -  Ij  +  I2  ♦  2Re<EL(t+T)E2(t)*>  (2.7) 

The  quantity  <E^ (t+t)E2(t)  >  is  defined  as  the  mutual  coherence  function 
ri2(t)  where  the  subscripts  denote  the  points  between  which  the 
coherence  is  measured. 

Complex  Degree  of  Coherence.  The  complex  degree  of  coherence, 
Y^2(t),  *8  3  normaH-ze<*  f°tm  of  the  mutual  coherence  function.  It  is 
defined  as 

Y12(t)  *  ri2<;T>/tril(0>r22(0>  1 1/2  (2,8) 

where  ru(0)  -  E1(t)E1(t)*  -  ^  and  ^(O)  -  E2(t)E2(t)*  -  I r  IQ 
can  now  be  written  as 

IQ  -  II  +  I2  +  21(1^2)  1/2  ]Re  y12(t)  (2.9) 

Visibility.  One  way  to  conduct  Young’s  two  slit  experiment  is 
under  what  Zernike  refers  to  as  best  conditions  (1:10).  These  are 
1^  -  I2  and  the  path  differences  are  small.  If  is  rewritten  as 

a  magnitude  times  a  phase,  i.e. 

Yi2(t)  “  t^r12(t)  !exp(i<ti12(T))  (2.10) 

where  <^2  represents  the  difference  in  phase  due  to  the  path  lengths, 
and  I...!  indicate  the  magnitude  of  the  quantity  they  enclose. 

Eq  (2.9)  may  be  rewritten  as 


IQ  -  211[1+Iy1-(t)!co8*i,(t)J 


(2.11) 


This  results  in  a  series  of  light  and  dark  fringes  appearing  on 
plane  B.  The  visibility  of  the  fringes,  V,  is  defined  as  (1:8) 

V  -  (I  -  I  .  )/(I  +  I  .  )  (2.12) 

max  min  max  min 

Applying  Eq  (2.12)  to  Eq  (2.11)  yields 

V  -!y12(t)1  (2.13) 

The  significance  of  this  result  is  that  the  modulus  of  the  degree  of 
coherence,  can  be  directly  related  to  the  measured  visibility 

of  the  fringes  (2:511).  This  result  will  prove  useful  later  in  this 
development. 

Quasl-monochromatlc  Light  Sources.  When  the  light  source  emits 
quasi-monochromatic  light,  r^^)  is  called  the  mutual  intensity 
function  and  is  denoted  by  J12.  The  complex  degree  of  coherence  is 
still  known  as  such  but  is  now  denoted  by  The  source  in  the 

optical  system  being  modeled  is  actually  non  quasi-monochromatic.  A 
filtering  operation  (described  in  the  section  on  calculating  the 
mutual  intensity  function)  takes  place  which  effectively  separates  the 
incoming  radiation  into  a  series  of  quasi-monochromatic  sources. 

van  Cittert-Zernike  Theorem.  Consider  the  situation  depicted  in 


Figure  2.3  below. 
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Fig  2.3.  Example  Situation. 

The  van  Cittert-Zernike  theorem,  then,  states  that 

...the  complex  degree  of  coherence,  which  describes  the  correlation 
of  vibrations  at  a  fixed  point  S2  and  a  variable  point  (see 
Figure  2.3)  in  a  plane  illuminated  by  an  extended  quasi- 
monochromatic  primary  source,  is  equal  to  the  normalized  complex 
amplitude  at  the  corresponding  point  S1  in  a  certain  diffraction 
pattern,  centered  on  S2«  This  pattern  would  be  obtained  on 
placing  the  source  by  a  diffracting  aperture  of  the  same  size  and 
shape  as  the  source,  and  on  filling  it  with  a  spherical  wave 
converging  to  S„,  the  amplitude  distribution  over  the  wavefront  in 
the  aperture  being  proportional  to  the  intensity  distribution 
across  the  source  (2:510). 

Therefore,  the  source  or  object  intensity  distribution  may  be  found 
by  taking  the  inverse  transform  of  the  complex  degree  of  coherence. 

This  is  the  basic  underlying  principle  on  which  the  optical  system  to  be 
modeled  operates.  One  last  useful  observation  is  that  the  sources 
considered  here  are  real  and  the  Fourier  transform  of  such  a  source  is 
complex  symmetric  due  to  the  hermitian  nature  of  the  transform  (6:193). 
This  means  that  if  the  complex  degree  of  coherence  p^2  is  measured  at  a 
particular  spatial  frequency  f  ^  then  it  is  automatically  known  at  the 
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symmetric  spatial  frequency  -f 


It  is  simply  the  complex  conjugate  of 


xl* 

the  value  of  at  ij* 

Linking  of  Theory  to  the  Optical  System 

The  optical  system  to  be  modeled  behaves  somewhat  like  a  coherent 
imaging  system.  According  to  Goodman  (7:106-110),  if  the  object  is 
illuminated  by  coherent  light  the  impulse  responses  comprising  the  image 
must  be  added  on  a  complex  amplitude  basis.  Therefore,  a  coherent 
imaging  system  is  a  linear  system  with  respect  to  complex  amplitude 
7:107).  In  a  coherent  imaging  system,  the  image  is  the  convolution  of 
the  image  predicted  by  geometrical  optics  with  an  impulse  response 
determined  by  the  exit  pupil  of  the  system  (7:105).  This  is  denoted  by 

00 

WV  *  ff  h(xi-xo*yi"yo)uo(xo’yo)dxodyo  (2*14) 

_oo 

where  u^  and  uq  are  the  image  and  object  amplitude  distributions  and 

00 

Mx^y^  -  ff  P(Xdix,A«|i1y)exp(-j2Tf(xix  +  yty))dxdy  (2.15) 

—oo 

and  where  P  is  the  pupil  function  and  d^  the  image  distance  (7:105).  Eq 
(2.15)  is  in  the  form  of  a  Fourier  transform  so  that  h  is  the  Fourier 
transform  of  the  pupil  function  P. 

Applying  the  convolution  theorem  of  Fourier  transforms  to  Eq  (2.14) 
yields 


G.(f  ,f  )  =  H(f  ,f  )G  (f  ,f  ) 
i  x  y  x  y  o  x  y 


(2.16) 


where  (denoting  a  Fourier  transform  by  F) 
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G^(f?,fy) 
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Eq  (2.16)  describes  a  linear  system.  Since  h(x,y)  is  the  Fourier 
transform  of  P(Ad^x.-Xd^y) ,  the  value  of  H(fx,f^)  is  P(-Ad^x,-Ad^y)  or, 
if  one  assumes  a  reflected  coordinate  system,  H(fx,f^)  ■  P(Ad^x,Ad^y) 
(7:110-111).  The  nature  of  H  is  that  it  allows  all  of  the  light  at  the 
sampled  frequencies  to  pass  through  and  completely  attenuates  the  light 
at  all  other  points.  Inverse  transforming  yields  the  image  u^. 

The  input  or  object  intensity  distribution  for  the  optical  system 

being  modeled  is  the  mutual  intensity  function.  H  consists  of  points 

sampled  by  the  optical  elements  of  the  system.  is  then  a  spatially 

filtered  version  of  Gq.  Inverse  transforming  G^  will  then  produce  the 

image.  H  determines  the  amount  of  spatial  filtering.  Consider  this 

one-dimensional  case.  If  G  was  a  rect(f  )  and  H  was  a  rect(2f  ) ,  G . 

o  x  x  1 

would  be  a  rect(2fx).  Spatial  frequencies  greater  that  0.5  are  filtered 
out.  This  obviously  degrades  the  quality  of  the  final  image.  This  is 
why  the  effect  of  the  aperture  or  pupil  function  is  of  such  Interest. 

Calculation  of  the  Mutual  Intensity  Function.  The  scene  of  in¬ 
terest  is  considered  to  be  spatially  incoherent,  temporally  stationary, 
and  non  quasi-monochromatic.  The  source  may  therefore  be  thought  of  as  a 
collection  of  m  independent  oscillators  operating  at  their  own  in¬ 
dividual  frequencies  and  radiating  a  complex  field  Em(t).  The  total 
field  present  at  each  of  lenses  and  L 2  can  then  be  thought  of  as  the 
sum  of  the  fields  due  to  each  oscillator.  This  is  denoted  by 

E^t)  -  EEml(t)  (2.17) 

m 

E2(t)  -  ^(t)  (2.18) 

m 


A  linear  filtering  operation  now  takes  place  in  order  to  separate 
the  Incoming  fields  for  detection  by  the  array  of  n  linear  detectors. 


The  filters  are  assumed  to  be  narrow  band  (to  fulfill  the  requirement  of 

quasi-monochromatic  light),  are  identical  for  each  lens,  and  have  a 

Fourier  transform  G  .  The  received  signal  is  then  Sf  (t)  -  E  (t)  * 
n  n  n 

g  (f  )  (n  ■  1,2,...)  (n  indicates  which  detector  is  being  analyzed), 
n  n 

The  mutual  intensity  function  can  then  be  rewritten  as  just 

r!2  '  ‘s,.l<t>st^«t>*>  '  Sfnl<t)Sfn2(t>'  <2‘19> 

This  illustrates  that  the  entire  continuous  mutual  intensity 
function  is  not  calculated  but  instead  only  at  the  n  specific  fre¬ 
quencies.  A  sampled  version  is  obtained  instead.  This  results  in  some 
degradation  of  the  image.  The  effects  of  the  number  of  frequencies 
sampled  on  image  quality  is  therefore  one  of  the  goals  of  this  thesis. 

The  Aperture  Function.  The  linear  detector  array  cannot 
(unfortunately)  contain  enough  detectors  to  detect  every  frequency 
present.  Therefore,  assume  that  the  detector  array  contains  four 
detectors  at  frequencies  f ^ ,  f2>  f^  for  the  following  example. 

Figure  2.1  illustrated  the  system  geometry.  The  mutual  intensity 
function  is  measured  at  intervals  along  £0  as  the  system  passes  across 
the  scene.  The  magnitude  of  the  four  frequencies  measured  by  the 
detector  array  is  recorded  for  each  interval.  A  polar  plot  of  the 
frequency  versus  A0  is  shown  in  Figure  2.4. 

The  system  is  able  to  obtain  resolution  in  the  dimensions  of  slant 
range  and  azimuth.  The  are  denoted  by  f ^  and  f^  in  Fourier  transforms 
as  in  Figure  2.4.  The  directions  of  the  radio  lines  can  be  thought  of 
two  ways  the  direction  is  mathematically  the  result  of  the  difference  of 
two  unit  vectors  denoting  the  positions  of  the  lens  relatively  to  the 
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Fig  2.4.  Polar  Plot  of  Frequencies  Sampled, 
source  (see  Appendix  A).  The  direction  can  be  thought  of  conceptually 
as  the  direction  of  the  lens  separation  relative  to  the  source  as  the 
lens  move. 

The  resolution  along  the  two  dimensions  is  determined  by  A0  the 

frequency  spacing.  This  is  best  illustrated  by  the  following  special 

cases.  When  the  lenses  are  (hypothetically)  infinitely  far  away  from 

the  source  such  that  they  lie  in  the  plane  of  the  source,  the  sampled 

frequencies  all  lie  along  f^  yielding  resolution  in  slant  range  but  none 

in  azimuth.  When  the  lenses  are  in  a  broadside  position  or  directly 

overhead,  no  resolution  along  ff  is  possible  because  all  the  sampled 

frequencies  lie  along  f  .  The  resolution  along  f  is  determined  by  the 

8  8 

spacing  between  samples. 

It  is  clear  from  Figure  2.4  that  not  all  frequencies  are  sampled. 
The  hermitian  nature  of  the  mutual  intensity  function  described  earlier 
allows  the  determination  of  the  function  p^  at  t*,e  points  -f^,  -^2’ 

-f^,  -f  ^  because  their  value  is  equal  to  the  complex  conjugate  of  p^  at 
the  points  f^,  f^»  f^.  However,  the  amount  of  frequency  coverage  is 
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limited  by  A0  and  the  number  of  detectors.  The  extent  of  A0  determines 
the  overall  shape  of  the  aperture  function  and  the  frequencies  at  which 
detectors  are  present  determine  which  frequencies  are  sampled  and  which 
are  not. 

Final  Output  of  the  Optical  System.  The  final  output  is  an  image 
of  the  original  source  obtained  by  inverse  Fourier  transforming  a 
sampled  version  of  the  mutual  intensity  function  as  governed  by  Eq 
(2.16).  This  results  in  the  degradation  of  the  image  since  not 
all  of  the  frequency  components  of  the  mutual  intensity  function  are 
present.  This  thesis  determined  what  effect  various  aperture  functions 
had  on  the  final  image. 
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III.  THE  COMPUTER  MODEL 


This  chapter  di-cusses  the  development  and  operation  of  the  com¬ 
puter  model  of  the  optical  system.  The  supporting  hardware  and  software 
will  be  described.  A  brief  explanation  of  the  overall  flow  of  infor¬ 
mation  will  be  given.  Appendix  B  contains  the  program  listings  and 
operating  instructions.  A  printout  of  a  sample  run  is  also  provided. 

Computer  Support 

The  model  was  developed  at  the  Electro-Optics  Branch  of  the  Air 
Force  Wright  Aeronautical  Laboratory.  The  computer  employed  was  a 
Hewlett-Packard  (HP)  9000  which  ran  under  the  UNIX  operating  system. 

The  graphics  support  consisted  of  an  HP  2623a  graphics  terminal  with  an 
internal  printer  and  an  HP  7550a  plotter.  The  computer  language  used  to 
construct  the  model  was  Fortran  9000  (a  version  of  Fortran  77).  The 
HP  9000  had  several  important  qualities  which  are  enumerated  below. 

One  of  the  best  features  was  that  of  being  a  virtual  memory 
machine.  This  eliminated  program  size  considerations.  The  only  concern 
was  speed  of  execution  since  the  machine  can  accomodate  a  program  of  any 
size.  This  was  one  less  problem  to  have  to  consider  and  therefore 
allowed  more  concentration  on  the  physics  of  the  thesis.  The  two 
megabytes  of  memory  were  more  than  sufficient  for  the  model. 

The  graphics  package  (named  Advanced  Graphics  Package  or  AGP)  was 
also  easy  to  use  (if  not  a  little  confusing  to  learn).  The  only  bad 
feature  was  that  there  was  no  way  of  hiding  lines  if  the  programmer  did 
not  know  what  kind  of  data  was  coming.  This  is  why  the  plots  that 
appear  later  in  this  thesis  are  only  one  quarter  of  the  front  half  of 
the  picture.  This  eliminated  many  of  the  confusing  lines  and  was 
possible  due  to  the  symmetry  of  the  information. 


The  feature  that  made  AGP  so  useful  was  the  .concept  of  a  Work 
Station  Program  (WSP).  Each  graphics  device  has  an  associated  WSP 


which  takes  care  of  all  device  dependent  affairs  allowing  the  programmer 
to  use  the  same  plotting  programs  on  different  devices.  The  only 
changes  occur  In  device  initialization  and  in  the  calling  sequences. 
Appendix  B  contains  more  details  on  this. 

The  Computer  Program 

The  preceding  flowchart  illustrates  the  various  parts  of  the 
computer  model  and  how  they  interact.  The  model  is  invoked  by  typing 
"synapt"  (SYNthetic  APerTure)  on  the  terminal.  This  is  also  the  name  of 
the  main  program.  See  Appendix  B  for  more  detailed  running  instruc¬ 
tions. 

The  model  first  initializes  all  graphics  devices.  It  is  then  ready 
to  find  out  what  type  of  source  is  to  be  Imaged.  The  model  contains 
five  preprogrammed  sources  and  will  allow  the  user  to  put  in  his  own. 

The  five  preprogrammed  sources  are  a  point  source,  a  two-point  source, 
an  edge,  a  slit,  and  a  circle.  The  model  analytically  computes  the 
Fourier  transform  of  the  first  four  sources  and  invokes  a  Fast  Fourier 
Transform  (FFT)  subroutine  to  compute  the  Fourier  transform  of  the 
circle  and  that  of  the  user's  own. 

The  model  then  asks  for  information  regarding  the  aperture  or  pupil 
function.  The  first  variables  required  are  the  range  to  the  source  and 
the  lens  separation.  These  variables  determine  what  frequencies  will  be 
sampled  (see  appendix  A  for  a  rigorous  derivation  of  this).  A  maximum 
A0  is  also  calculated  from  the  range  information,  collector  speed,  and 


collector  stability  information  provided  by  the  user.  The  user  then  is 
asked  for  A©.  An  upper  and  lower  limit  are  required.  This  A©  is 
compared  to  the  maximum  possible  AS  calculated  above.  If  the  calculated 
AS  is  exceeded,  the  user  must  start  over  with  new  range  and  lens 
separation  information. 

It  was  decided  to  use  rectangular  symmetry  instead  of  radial 
symmetry  in  the  aperture  and  source  distribution  because  an  FFT  sub¬ 
routine  that  worked  with  radial  symmetry  could  not  be  found.  The  points 

in  the  pupil  function  were  modelled  as  having  a  radial  distance 

2  2  1/2 

(  (f  +  f  )  where  ff  and  ffi  are  the  rectangular  components  of  the 
frequency  in  range  and  azimuth  respectively)  which  was  the  equivalent 
to  the  radial  frequency  at  that  point.  Each  element  of  the  pupil 
function  either  transmits  entirely  or  attenuates  entirely.  A  picture  of 
the  pupil  function  appears  in  Figure  3.2. 

All  of  these  frequencies  are  not  sampled.  The  frequencies  that  are 
sampled  are 

f  -  2(f  /c)  sin  (a/2)  (3.1) 

m 

where  f^  is  the  frequency  of  a  filter  and  a  is  the  angle  between  the 
lenses  formed  by  and  R£  (see  appendix  A  for  supporting  information 
and  a  derivation  of  Eq  3.3).  The  spatial  frequencies  f  that  are  sampled 
in  the  optical  system  being  modelled  are  0  to  64  1/m  based  on  ranges  of 
1  to  3  kilometers,  a  lens  separation  of  0.5  meters,  a  wavelength  (c/f^) 
band  of  8  to  12  pmeters,  and  a  maximum  pupil  radius  of  16  bits  in  a  256 
by  256  array  of  points  that  modelled  the  pupil.  The  model  computes  the 


upper  and  lower  frequencies  based  on  the  information  Indicated  above  and 
shades  the  appropriate  area  of  the  pupil  to  reflect  the  operating 
conditions. 

The  source  FFT  now  undergoes  an  inverting  process.  The  FFT  sub¬ 
routine  normally  places  the  high  frequency  components  in  the  middle  of 
the  transformed  array  and  the  low  frequency  components  at  the  corners. 
This  process  allows  the  FFT  to  appear  in  its  more  commonly  recognized 
form  with  the  low  frequencies  in  the  center  and  the  high  frequencies  on 
the  outside  edges.  See  Appendix  B  to  see  how  it  is  implemented. 

The  source  FFT  now  is  multiplied  by  the  pupil  function.  The  numbers 
to  be  multiplied  in  general  are  complex  and  the  operation  is  of  the  form 

(a  +  ib)  (  c  +  id  )  *=  Answer  (3.2) 

where  i  is  the  square  root  of  negative  one,  a+ib  represents  the  value 

of  a  point  of  the  FFT  of  the  source,  and  c  +  id  is  the  value  of  a  point 

in  the  pupil.  However,  since  the  components  of  the  pupil  function  are 
entirely  real  (  d  -  0  ) ,  the  result  is 

ac  +  icb  =  Answer  (3.3) 

This  result  is  stored  in  the  array  that  originally  held  the  source  FFT. 

The  final  answer  is  now  obtained  by  taking  the  inverse  Fourier 

transform  (IFT)  of  the  result  of  Eq  3.3.  The  data  is  inverted  as  before 
in  order  for  the  answer  to  appear  in  its  original  form.  The  user  can 
get  a  picture  of  this  and  compare  it  to  the  original  source  to  see  just 
how  well  his  system  has  performed. 

There  are  two  subroutines  other  than  the  main  program  which  do  a 
majority  of  the  work.  The  most  important  is  the 


FFT  subroutine.  Figure  3.3  compares  the  results  of  an  analytically 


computed  32  sine  (32  f^)  to  the  Fourier  transform  of  a  pulse  with  a 
width  of  32. 

The  subroutine  computes  the  transform  by  finding  the  terms  of  a 
Fourier  series  of  the  same  function  as  if  it  were  indeed  periodic.  A 
period  of  236  units  was  empirically  found  to  yield  an  acceptable  ac¬ 
curacy  as  illustrated  in  Figure  3.3.  Although  a  larger  period  was  more 
desirable,  lengthening  the  array  also  resulted  in  greatly  Increasing  the 
time  required  by  the  computer  to  compute  the  FFT.  The  subroutine 
required  the  input  of  both  a  real  and  imaginary  component.  This  re¬ 
sulted  in  requiring  two  256  by  256  element  arrays  to  adequately  describe 
the  source  in  two  dimensions.  A  listing  of  the  subroutine  may  be  found 
in  Appendix  B. 

The  other  subroutine  is  called  PLTRST  for  PLoTReSulT.  This  sub¬ 
routine  carries  out  the  graphics  operations.  The  subroutine  requires 
that  the  calling  program  indicate  which  plot  is  needed.  The  subroutine 
then  draws  the  appropriate  axes  in  the  appropriate  projection  and  types 
all  of  the  appropriate  headings,  titles,  and  other  markings.  This  is 
output  to  the  appropriate  device  as  the  user  has  indicated.  The  WSP 
takes  care  of  the  actual  drawing.  See  Appendix  B  for  more  details. 

Several  other  subroutines  are  also  invoked  in  the  model.  However, 
they  perform  only  support  functions  and  it  is  therefore  unnecessary  to 
go  into  detail  here  on  what  they  do.  Again,  Appendix  B  contains  more 
information  on  this  matter. 

Outputs 

The  outputs  of  the  model  are  all  graphical  in  nature.  The  model 
will  plot  on  both  the  terminal  and  the  plotter.  There  are  a  total  of 


five  possible  outputs  on  any  one  run  of  the  model.  These  are  (1)  of  the 
source  irradiance  distribution,  (2)  the  aperture  transmittance  function, 
(3)  the  FFT  of  the  source,  (4)  the  product  of  the  source  and  the  aper¬ 
ture,  and  (5)  the  inverse  FFT  of  the  product  of  the  source  and  the 
aperture  or  the  image. 

The  scales  in  the  object  and  the  image  plane  are  in  terms  of  a 
dimensionless  variable  V  defined  as 

V  **  (2xa/(Xdj))x^  (3.4) 

where  a  is  the  pupil  radius,  X  is  the  cutoff  wavelength,  d^  is  the  image 
distance,  and  is  the  position  of  a  point  in  the  image.  The  size  of  a 
is  0.5  m  which  was  found  from  the  relationship  found  in  reference  7:112 
for  the  cutoff  frequency  of  a  coherent  transfer  function 

f  -  1/ (2Xd.)  (3.5) 

c  1 

where  1  is  the  diameter  of  the  pupil.  Equating  this  to  Eq  3.1  resulted 
in  1  -  d  where  d  is  the  lens  separation.  The  value  of  d^  is  fixed  due 
to  the  fixed  focal  length  of  the  lenses.  The  value  of  X  is  taken  to  be 
8um  at  the  radius  of  the  pupil.  A  sample  object  size  of  3.125  cm  was 
used  to  arrive  at  a  corresponding  value  of  V  of  2000n/512  .  This 
object  wa6  assumed  to  be  32  bits  wide  in  the  256  by  256  array  represent 
ing  the  source.  As  noted  earlier,  only  one  quarter  of  the  front  half  of 
the  plot  is  illustrated.  Therefore,  the  value  of  V  must  be  doubled  to 
find  the  total  width  or  length  of  the  object.  The  tick  marks  in  the 
plots  of  the  object  and  the  image  plane  represent  8  bits  in  the  source 
array.  This  is  why  the  scale  appears  as  is  illustrated  in  Fig  3.4. 
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The  other  outputs  occur  in  the  frequency  plane.  The  scale  is  as 
described  above  in  the  pupil  function.  The  axes  are  labelled  in  range 
(f^)  and  azimuth  (f^).  A  sample  of  the  modulus  versus  the  spatial 
frequencies  of  the  FFT  of  a  6.25  cm  wide  edge  appears  in  Fig  3.5. 

All  of  the  plots  are  normalized  in  amplitude.  In  all  cases,  the 
modulus  of  the  amplitude  is  plotted.  The  factor  by  which  the  informa¬ 
tion  is  normalized  is  shown  on  the  plots  as  Rnorm.  The  preprogrammed 
sources  are  assumed  to  have  an  amplitude  of  one  originally. 


IV.  RESULTS 


This  chapter  contains  the  results  of  successive  runs  of  the  com¬ 
puter  model.  Five  types  of  sources  were  imaged  through  five  types  of 
apertures.  The  results  are  in  the  form  of  plots  of  the  image  of  each  of 
the  five  sources  due  to  each  aperture  configuration.  Each  aperture  and 
the  images  of  the  five  sources  due  to  each  aperture  appear  on  foldouts 
at  the  end  of  the  chapter.  R.  Barakat  (16:205-223)  has  also  examined 
the  effects  of  different  apertures  on  the  images  of  various  objects. 

His  results  agree  quite  closely  with  the  results  obtained  with  the 
computer  model. 

The  initial  aperture  reflected  the  best  case  scenario;  i.e.  all 
frequencies  were  allowed  to  pass  below  the  cutoff  wavelength  of  8  pm. 

The  next  two  apertures  reflect  how  the  images  degrade  as  fewer  and  fewer 
of  the  low  frequencies  were  allowed  to  pass.  The  images  began  to  show 
the  characteristics  of  edge  enhancement  (10:61-62)  as  the  lower  fre¬ 
quencies  were  filtered  out.  These  apertures  established  a  baseline  from 
which  apertures  obtained  under  realistic  conditions  could  be  compared. 

Figures  4.1  and  4.2  illustrate  the  process  of  edge  enhancement  for 
the  case  of  an  edge  source  with  a  uniform  amplitude  distribution  and  a 
width  of  6.25  cm.  The  X's  in  Figure  4.1  indicate  the  bounds  on  the 
spatial  frequencies  of  the  Fourier  transform  of  the  edge  which  are 
allowed  to  pass  by  the  actual  system.  This  figure  reflects  a  lens 
separation  of  0.5  m  and  a  range  of  1  km.  These  two  factors  plus  the  8 
to  12  pm  bound  on  the  detectable  wavelengths  combine  to  yield  a  limit  on 
the  spatial  frequencies  of  approximately  40  to  64  cycles  per  meter. 

Figure  4.2  is  the  image  obtained  by  analytical  methods  with  the 
above  mentioned  system  parameters.  The  edge  falls  at  X^  equal  4  in 

Figure  4.2.  This  confirms  that  edge  enhancement  takes  place. 
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Figures  4.3  and  4.4  compare  how  the  computer  model  behaves  under 
identical  operating  conditions.  Figure  4.3  illustrates  which  fre¬ 
quencies  are  sampled  by  the  computer  model.  The  analytical  case  above 
assumed  all  spatial  frequencies  between  40  and  64  cycles  per  meter  were 
sampled.  Figure  4.3  shows  that  this  is  not  the  case  in  the  computer 
model.  The  bounds  on  the  spatial  frequencies  remained  the  same  but 
samples  were  taken  only  at  4  cycles  per  meter  intervals  as  Indicated  by 
the  X's  in  Figure  4.3.  Figure  4.4  shows  the  image  obtained  by  the 
passive  system  through  the  computer  model  when  only  discrete  samples 
of  the  Fourier  transform  are  taken.  Figure  4.4  agrees  very  closely  with 
Figure  4.2  which  would  be  the  best  image  that  could  be  obtained.  This 
comparison  was  the  final  test  in  validating  the  computer  model. 

The  five  sources  imaged  were:  a  point  source,  two  point  sources 

with  a  separation  in  of  4,  an  edge  (that  was  modeled  by  a  rect  with  a 

width  of  4  along  V  ),  a  slit  with  dimensions  of  (V  X  V  )  4  X  2,  and  a 
fl  r  d 

circle  with  a  radius  of  4.  The  amplitude  distribution  was  unity  at  all 
points  on  the  sources.  Figures  4.5  and  4.6  illustrate  these  objects. 

The  Fourier  transforms  of  these  sources  were  taken  as  described  in 
Chapter  III  and  passed  through  five  different  apertures.  The  resultant 
images  of  these  sources  through  each  of  the  apertures  appear  on  five 
foldouts  at  the  end  of  this  chapter  along  with  the  aperture  used.  The 
five  cases  are  described  below. 

Case  I 

The  first  case  considered  what  images  could  be  expected  from  a  full 
aperture.  The  highest  frequency  present  was  64  cycles  per  meter.  The 
images  are  as  expected  with  some  ringing  present.  The  impulse  response 


A.  *.  *. 
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is  good  with  low  sidelobes.  The  images  can  be  interpreted  easily  at 
this  point.  This  case  is  reflected  on  the  first  foldout  (Figure  4.7). 


Case  II 

The  second  case  used  an  aperture  that  passed  only  spatial  fre¬ 
quencies  between  16  and  64  cycles  per  meter.  The  impulse  response  now 
is  more  spread  out  with  higher  sidelobes.  The  images  are  no  longer 
easily  interpreted.  Substantial  ringing  is  beginning  to  occur.  This 
case  is  illustrated  on  the  second  foldout  (Figure  4.8). 

Case  III 

This  case  reflects  what  happens  to  the  images  when  only  the  spatial 
frequencies  between  40  and  64  cycles  per  meter  are  sampled.  This  case 
reflects  realistic  spatial  frequencies  since  they  reflect  realistic 
system  parameters  of  0.5  m  for  lens  separation  and  1  km  for  range.  A0 
is  considered  to  be  180  degrees  in  this  case  still.  The  images  have 
degraded  even  more  with  more  ringing  present  along  V  .  This  is  re- 
fleeted  in  the  third  foldout  (Figure  4.9). 

Case  IV 

This  case  reflects  the  system  performance  under  the  conditions  in 
Case  III  with  the  addition  of  a  realistic  A 9.  This  case  is  a  reflection 
of  the  kind  of  data  that  can  be  expected  from  a  passive  synthetic 
aperture  system.  The  velocity  and  stability  used  in  calculating  A0  were 
880  ft/s  and  10  s.  This  resulted  in  a  A0  of  approximately  106  degrees 
which  was  centered  on  both  sides  of  fg.  The  impulse  response  has  even 
higher  sidelobes  and  the  images  have  become  slightly  noisier  than  in  the 
previous  case.  The  data  for  this  case  is  shown  in  the  fourth  foldout 
(Figure  4.10) . 
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The  last  case  was  Intended  to  show  what  happens  under  a  slightly 
different  system  configuration.  The  lens  separation  has  been  halved  to 
0.25  m.  The  sampled  spatial  frequency  range  is  now  approximately  20  to 
32  cycles  per  meter.  The  impulse  response  now  exhibits  wider  and  more 
spread  out  sidelobes  along  Va  than  in  previous  cases.  This  is  also 
reflected  in  the  images  of  the  other  sources.  This  is  illustrated  in 
the  last  handout  (Figure  4.11). 
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V.  CONCLUSIONS 


The  chapter  summarizes  the  conclusions  reached  on  the  performance 
of  the  passive  synthetic  aperture  system  and  the  computer  model.  These 
conclusions  are  based  on  the  results  presented  in  the  previous  chapter 
and  on  the  vast  amount  of  experience  gained  in  the  actual  operation  of 
the  computer  model. 

Conclusions  on  Optical  System 

1 .  The  system  will  behave  like  a  high  pass  filter.  This  is  because 
the  DC  component  of  spatial  frequency  will  never  be  measured  under 
realistic  operating  conditions. 

2.  The  images  will  be  edge  enhanced.  This  is  due  to  the  high  pass 
nature  of  the  system.  The  images  will  not  simply  be  a  slightly  degraded 
version  of  the  geometric  image  but  will  be  highly  complex  instead. 

3.  A  large  AO  is  desirable  to  improve  resolution  along  Va.  Limiting 
0  limits  resolution  along  Va  accordingly. 

4.  New  methods  of  image  interpretation  will  need  to  be  developed  in 
order  for  the  system  to  be  usable.  The  reasons  cited  in  conclusion  two 
will  require  the  development  of  new  techniques  and  algorithms  in  order 
to  interpret  the  information  correctly.  This  should  be  possible  since 
the  impulse  response  for  any  system  configuration  or  change  in  operating 
conditions  can  be  found.  Knowing  the  impulse  response  of  the  system 
should  enable  one  to  find  the  image  resulting  from  any  input. 

Conclusions  on  Computer  Model 

1.  The  model  performs  satisfactorily.  The  results  have  been  ver¬ 


ified  analytically  as  correct. 


2.  A  radial  FFT  program  would  allow  for  the  exact  simulation  of  the 
system.  The  sampled  spatial  frequencies  are  actually  distributed  In  a 
radial  fashion  (see  Appendix  A) .  Fitting  a  rectangular  grid  to  the 
situation  was  an  Interim  and  time  saving  solution.  A  more  exact  sim¬ 
ulation  could  be  obtained  through  the  use  of  radial  coordinates.  A 
program  that  takes  an  FFT  in  radial  coordinates  could  not  be  found. 
Therefore,  the  rectangular  grid  was  fitted  in  order  to  utilize  con¬ 
ventional  FFT  programs. 

3 .  Accuracy  can  be  increased  by  increasing  the  array  size.  This 
will  result  in  more  samples  of  the  input  and  its  Fourier  transform 
resulting  in  a  more  exact  representation.  This  will  also  increase 
program  size  and  slow  down  processing  considerably. 

4.  Hiding  lines  on  plots  would  increase  the  usefulness  of  the  data. 
However,  the  method  for  doing  this  is  not  readily  apparent. 

5.  The  effects  of  phase  should  be  determined.  A  phase  term  is 
present  in  the  mutual  coherence  function  as  shown  in  Appendix  A.  Its 
effects  have  not  been  studied  in  this  thesis. 
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Appendix  A 

This  appendix  relates  the  simple  case  of  a  scene  rotating  beneath 
the  lens  system  to  the  more  complicated  case  of  the  lens  system  being 
carried  on  a  collection  platform  and  moving  past  the  scene.  This  will 
require  a  rigorous  derivation  of  the  propagation  of  the  mutual  coherence 
function.  The  following  paragraphs  will  present  needed  background  and 
terminology  for  the  derivation  which  follows.  The  material  presented  in 
this  chapter  is  extracted  mainly  from  reference  14:2-4  to  2-13.  The 
system  geometry  is  illustrated  again  in  Figure  A.  1  below. 
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Fig  A.l.  System  Geometry. 


The  major  assumptions  are  that  the  measurement  (slant  range)  plane 


is  essentially  the  same  as  the  ground  plane  being  imaged.  The  lens 
separation  is  constant  and  the  system  moves  in  a  direction  parallel  to 
this  separation. 


The  first  step  is  to  derive  the  mutual  coherence  function  at  the  source 
and  next  to  propagate  it  to  the  lens  system. 

The  source  to  be  imaged  has  a  field  amplitude  E(r_'  ,t)  which  is 
spatially  incoherent  and  temporally  stationary.  The  field  at  any  point 
is  uncorrelated  with  any  other  point.  Therefore,  the  mutual  coherence 
function  as  evaluated  at  two  points  on  the  source  is 

ri2<T)  “  <E(£i'»t1>E(E2,,t2)*  >  (A,1) 

Ti2 (t)  can  be  reduced  as  in  Eq  2.6.  This  yields 

T12(t)  =  <E(r1,,t)E(r2',t-T)*  >  (A.2) 

Rewriting  Eq  A.2  in  terms  of  intensity  yields 

T12(t)  -  I<r,,T)«(ri»  -  r2')  (A.3) 

where  6  is  a  dirac  delta  function  denoting  the  spatial  incoherence  of 
the  source. 

Refer  now  to  Figure  A.2. 


Fig  A.2.  Vector  relationships. 


The  magnitude  of  the  vectors  r.  and  r -  denote  the  distance  from  the 

1  A  1  A 


center  of  the  source  to  each  lens,  r^  and  r^  are  unit  vectors  along  r^ 
and  r^,  and  jr_*  is  the  position  vector  of  a  point  in  the  scene  to  be 
imaged.  The  field  present  at  the  lenses  can  be  found  through  Huygens 


principle.  This  is  (within  a  constant) 


E(r,t)  -  /  E(r/,t  -  R/c)  drj /R 
A 


(A.  A) 


where  R  ■  ! _r  -  r_* ! ,  c  is  the  speed  of  light,  and  R/c  is  the 


lag  time  from  r_f  to  £. 


This  received  signal  (s  )  at  each  lens  passes  through  a  linear 

m 


filter  (separate  but  identical  for  each  lens)  with  a  center  frequency 


f^.  This  process  can  be  denoted  by 


s  (r,t)  -  E(r,t)  *  h  (t) 
m  —  —  m 


(A.5) 


where  hm(t)  is  the  electronic  impulse  response  of  the  filter  and  * 


denotes  a  convolution.  The  subscript  ra  denotes  the  mth  filter  in¬ 


dicating  that  m  total  frequencies  are  being  sampled  at  a  given  time. 


A  rigorous  derivation  of  the  propagation  of  the  mutual  coherence 


function  is  now  made.  The  theoretical  basis  may  be  found  in  reference 


2:537-599.  The  mutual  coherence  function  is  evaluated  in  terms  of  the 


received  signal  at  each  lens.  This  is  written  as 


r (— i  »  i2  *  T)  "  <8ml(— l' ,tl)  Sm2 (— 2 ' ,t2)  > 


(A. 6) 


where  sm^(r^,t^)  is  the  received  signal  at  lens  i  at  the  mth  frequency 


received  at  a  time  t^.  can  be  rewritten  by  substituting  for 


'v-o: 'a 


.*-1 
ii 


-zm 


mi 
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8  .  (r . ,t.)  as  allowed  by  the  relationship  in  Eq  A. 5.  This  results  in 
ml  1  1 


T12(t)  -  <(E(r1'.t1-R1/c)*hm(t1))(E(r2',  t2-R2/c)*hm(t2))  >  (A.7) 


This  equation  can  be  rewritten  once  more  making  a  substitution  for 


E(ri,ti)  as  allowed  by  Eq  A. A.  This  results  in 


T12(t)  =  </((  E(r1',t1-R1/c)dr,/Rl)  *  h^)) 

A 


/((E(r2’,t2-R2/c)dr'/R2)  *  hm(t2))  >  (A.8) 

A 


Replacing  the  convolution  symbol  with  the  convolution  integral  and 
removing  the  integrals  in  jr'  and  the  constants  R^  and  R2  outside  the 
time  average  (since  they  have  no  time  dependence)  yields 


ri2(x)  -=  1/(R1R2)  //</(E(r1,,t,-R1/c)hm(t1-t,))df)/(E(r2’,t"-R  /c) 
AA  -®  *  ^ 

h  (t0-t")dt")  >dr 'dr '  (A.9) 

m  A  —  — 


where  t'  and  t"  are  dummy  variables  introduced  by  the  use  of  the  con¬ 
volution  integral. 

Since  the  time  average  only  applies  to  E  because  of  its  rapid 
oscillations „  Eq  A.9  can  be  rewritten  as 


T12(t)  -  1/(R ^2)HIS  <  (E(r1',t'-R1/c)E(r2',t"-R2/c)  >  hjt^t') 
AA-®  . 


h  (t~-t")  dt'dt"  dr ' dr '  (A. 10) 

m  z  —  — 


Employing  Eq  A. 3  allows  the  quantity  within  the  time  average  brackets 
<  and  >  to  be  rewritten  which  yields 


T12(t)  -  l/(R1R2)////I(r',t'-t"-(RrR2)/c)6(r1*-r2,)dr'  h  (tj-t*) 
AA-®  ^ 

h  (t„-t")  dt'dt"dr'  (A. 11) 


The  delta  function  results  In  the  elimination  of  one  area  integral  which 
reduces  Eq  A. 11  to 


ri2(t)  -  l/(R1R2)///I(r',t'-t"-(R1-R2)/c)hm(t1-t,)hm(t2-t")  dt’ 
A-00 


*it"dr '  (A.  12) 


The  next  step  is  to  rewrite  the  filter  impulse  responses  as  their 
Fourier  transforms  which  are  denoted  Hm(f)  (the  same  symbol  is  used  for 
both  filters  since  they  are  identical).  These  are  (in  integral  form) 


h  (t.-t')  «/  H  (f)exp(i2xf(t,-t'))df 
ml  m  i 


(A. 13) 


h(t9-t")  -/  H  (f)  *  exp(-i2uf (t_-t"))df 

m  i  m  / 


(A. 14) 


The  product  of  Hm(f)  and  Hffl(f)  may  be  taken  outside  the  integrals  since 
they  are  constants  for  all  frequencies.  Replacing  the  filter  impulses 
in  Eq  A. 12  with  their  Fourier  transforms  of  Eqs  A. 13  and  A. 14  results  in 


ri2(t)“,Hm(f)  )1(R1R2)///  I(£,.t'-t"-(R1-R2>/c)/exP(127rf(ti’,:,>>df 

A— 00  00  -00 

/exp(-i2irf(t2-t”))df  dt'dt'Mr'  (A.  15) 

-00 

The  integrals  of  the  exponential  functions  are  in  the  forms  of  delta 
functions.  Therefore,  Eq  A. 15  can  be  rewritten  as 


F12(T)  “  >Hm(f)^!/(R1R2)///I(r',tf-t"-(R1-R2)/c)6(t1-t') 

A-00 

<5(t„-t")  dt’dt’Mr'  (A.  16) 


Carrying  out  the  integrations  with  respect  to  t'  and  t"  utilizing 
the  sifting  property  of  the  delta  function  yields  the  following  result 


ri2(t)  -  (!Hm(f)!i)/(R1R2)/l(r’,t1-t2-(R1-R2)!c)  dr*  (A. 17) 


Denoting  t^-t2  as  T  and  rewriting  I  as  a  Fourier  transform  results  in 
the  following  result  for  the  mutual  coherence  function 


r,2(T)  -  l/U^R  )//!H  (f)!2l(r',f)exp(i2TTf(T-(R1-R2)/c))dfdr'  (A.18) 


I(£*  ,f)  is  essentially  a  constant  with  respect  to  f  at  the  fre¬ 
quencies  of  concern.  I(_r' ,f)  will  therefore  be  denoted  as  Kir' )  from 

this  point  on.  It  is  also  assumed  that  H  (f)  is  sufficiently  narrowband 

m 

that  the  factor  of  f  in  the  exponent  in  Eq  A. 18  can  be  replaced  by  its 


center  frequency  value  f  .  The  actual  bandwidth  of  the  filter,  B,  is 

m 


defined  as 


B  *>/!H  (f ) !  df 
_»  m 


(A.  19) 


Lastly,  a  far  field  assumption  is  made;  i.e. 

A 


Rl,2  "  -1,2  "  rl,2  -1,2 


(A. 20) 


where  r.  _  is  the  magnitude  of  jr. 

1  9  /  1 


The  final  form  of  r^2(-t)  with  these  assumptions  is 


ri2 (T)  “  C/I(_r')exp(-i2xfm(r2-£1)  ,£,/c)d_r' 
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(A. 21) 


where  C  ■  (B/R,R„)exp(-i2nf  (T-(r_-r, ) /c) .  The  quantity  to  be  inte- 
I  l  m  —L  — i 

grated  in  Eq  A. 21  is  in  the  form  of  a  Fourier  transform.  This  is 
denoted  symbolically  as 


r  1 2 )  “  (B/R1R2)exp(-i2nfj|i(T-(ri2-r1)/c)*(f) 


(A. 22) 


where  f  ■  (f  /c) (r~  -  r,).  If  a  is  the  angular  separation  of  the  lenses 
—  m  — L  — l 

(and  therefore  the  angle  between  r^  and  x_^) ,  £_  can  be  rewritten  as 


f  **  (2f  /c)  sin  (a/2)  f 
—  m 


(A. 23) 


where  f  is  a  unit  vector  in  the  direction  of  r_  -  r, .  Vector  algebra 
A  z  l 

shows  that  f  is  also  perpendicular  to  the  bisector  between  and  r^ 
The  case  of  a  rotating  scene  may  now  be  considered.  Figure  A. 3 
illustrates  the  current  situation. 


A© 


All  variables  are  as  defined  In  earlier  chapters.  Note  the  unit 
A  A 

vector  f.  Any  frequency  along  this  unit  vector  (denoted  f.  ■  f^f  )  can 

be  measured  by  either  varying  the  center  frequency  f  or  by  placing  m 

linear  detectors  in  the  system.  The  m  linear  detectors  would  allow  the 

simultaneous  measurement  of  at  m  different  frequencies  as 

governed  by  Eq  A. 22.  The  measurements  are  made  at  specific  intervals  as 

the  scene  rotates  until  a  total  rotation  of  A0  ■  180  degrees  has  been 

made.  No  more  samples  need  be  taken  since  the  transform  being  sampled 

is  hermitian  (see  Chapter  II) . 

Now  consider  the  situation  in  Figure  A. 4  where  the  lens  system  is 
moving  while  the  scene  remains  stationary. 


The  situation  is  obviously  identical  to  the  rotating  scene 

A 

scenario.  The  movement  of  the  lens  system  allows  the  unit  vector  f 
to  sweep  out  over  an  angle  of  A0  *  180  degrees  just  as  in  the  case  of 
the  rotating  scene.  The  m  frequencies  can  be  measured  as  in  the 
stationary  scene  case.  Figure  A. 4  shows  where  the  samples  are  taken  as 
does  Figure  A. 3.  . The  frequencies  measured  for  the  four  values  of  A0 
shown  are  the  same  in  both  cases.  A  polar  plot  of  A0  versus  the  spatial 
frequencies  f  and  f  for  both  cases  yield  the  same  results. 


Fig  A. 5.  Polar  Plots  of  A0  vs  f  and  f  for  Rotating  and  Stationary 
Cases.  r  a 


This  appendix  contains  the  source  listings  of  the  computer  model 


programs,  and  a  sample  model  run  illustrating  how  to  use  the  model. 
Computer  Listings 

The  listings  presented  in  the  following  pages  are  in  the  order  in 
which  they  occur  when  the  model  is  run.  The  programs  are  written  in 
Hewlett-Packard  Fortran  9000  and  are  heavily  commented. 
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c  *  »  * 

1  0 

c 

1 1 
12 

13 

14 

c 

15 

C*  t  * 

16 

c 

17 

c 

18 

c 

19 

c 

2  0 

c 

21 

c 

22 

c 

23 

c 

24 

c 

25 

CM* 

26 

c 

27 

28 
29 

c 

30 

C  *  *  * 

31 

c 

32 

c 

33 

c 

34 

c 

35 

Cm  (■ 

36 

c 

37 

38 

39 

c 

40 

c  ♦  *  *  ' 

4  1 

c 

42 

c 

43 

c 

44 

C  *  *  *  i 

45 

c 

46 

10 

47 

c 

49 

C  MM  i 

49 

c 

50 

c 

51 

c 

52 

c 

53 

C  ♦  M  1 

program  s  y  nap  t 


This  is  the  main  program.  The  relevant  data  is  input  via 
prompts  by  the  program.  The  program  also  prompts  the 
operator  as  to  what  kinds  of  output  he  desires. 


Dimension  Real  <  256  )  ,  Rimag  <  256  )  ,  RSourc  <  256 
iCSourc  <  256  ,  256  )  ,  RApert‘<  256  ,  256  )  ,  Ar  (  2 
26  ,  256  ) 


256  )  , 
Plot  (  25 


These  Dimension  statements  declare  the  arrays  Real,  Rimag, 

<  which  are  used  to  pass  data  to  the  FFT  subroutine  ),  RSourc 
and  RApert  (  which  contain  the  real  components  of  the  source  and 
aperture  functions  respectively  )  and  CSourc  <  which  contains 
the  imaginary  part  of  the  source  distribution).  The  program 
prompts  the  operator  for  the  source  amplitude  distribution 
and  the  aperture  function  of  interest. 


Common  /  Args  /  Real 
1,  Plot  ,  Ar  ,  Lower 


Rimag 
I  upper 


RSourc 

Length 


CSourc 

Iwide 


RApert 


This  common  area  will  be  used  to  pass  arguments  more  ef¬ 
ficiently  to  the  subroutines  that  require  them. 


Character  Answer 
Data  y/'y'/ 


Prompt  operator  for  type  of  aperture  and  source  distributions. 


Call  Grphin 


Get  the  source  type.  Unit  7  is  the  screen  and  unit  5  the 
keyboard  . 


**tMMMMMMMMM*M*M*MM*MMM** 
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54  C 

55 

56  30 

57 

58 

59 

60 
6  1 
62 

63 

64 

65 
6  6 

67  C 

6  8  cm 


Write  <  7  ,  30  ) 

F  o  r  m  a  t  < /  ,  'Enter  your  source  irradiance  distribution.  You  nay  choo 
lse  fron  one  of  the  pre-  programmed  distributions  below  or  create  y 
2our  own.  Type  the  numiber  of  your  'selection  after  finding  you 

3r  choice  on  the  Menu  below  when  proMpted. 

4'A  point  source T35  1 ',/ , 

5 ' A  two  point  source :  '  ,  T  3  5  , '2',/. 

6 'An  edge:',T35,'3',/f 
7 ' A  slit :  '  ,  T35  ,  ' 4 '  ,/ , 

5'A  circle  of  variable  r ad i us : ' , T35 , ' 5 ' , / , 

6/Your  own  cr e a t i on T 35  6 ',/ , 

7/, %,  'Enter  your  selection  [1-61:  ') 


vs;*. 


< 


69 

70 

71 

72 

73 

74 


C 

C 

C 

C 

C 

c 


Call  the  appropriate  subroutine  to  set  the  appropriate  entries 
in  RSourc  and  CSourc  to  reflect  the  type  of  source  desired  if 
the  source  chosen  is  not  one  of  the  five  preprogrammed  ones. 


t ************* 


7  5  C 


76 

Read  ( 

5  ,  40 

) 

Ichose 

77 

4  0 

F  o  r  m  a  t 

(ID 

78 

If  < 

Ichose 

.  e  ( 

1 

) 

Goto 

45 

7? 

If  < 

Ichose 

•  eq 

2 

) 

Call 

T  w  o  p  n  t 

30 

If  < 

Ichose 

eq 

.  3 

) 

Call 

Edge 

31 

If  < 

Ichose 

eq 

4 

) 

Call 

Slit 

82 

If  < 

Ichose 

eq 

.  5 

) 

Call 

Circle 

33 

If  ( 

Ichose 

■  eq 

.  6 

) 

Call 

Other 

8  4  C 

85  C**M**»* *************************************************************** 

86  C 

87  C  See  if  operator  needs  a  picture  of  the  source  being  Modelled. 

88  C 

8?  C************* ****************************************  ****************** 

9  0  C 

91  45  Write  (  7  ,  50  ) 

92  50  Fornat (/ , 'Do  you  want  the  source  plotted  on  the  screen  I  y/n  ]  ? 

9  3  ID 


94  Read  (  5  ,  60  )  Answer 

?5  60  Fornat(Al) 

?  6  If  (  Answer  ,ne.  y  )  Goto  65 

?7  Call  Pltrst  <  1  ,  Ichose  ) 

?8  65  Write  (  7  ,  70  ) 

?9  70  Fornat(/,#,'Do  you  want  a  hardcopy  of  the  results  C  y/n  1  ?  ') 

00  Read  (  5  ,  60  )  Answer 

01  If  (  Answer  .ns.  y  )  Goto  75 

02  Call  Plot  in 

03  Call  Pltrst  <  1  ,  Ichose  ) 

04  Call  Plotof 

0  5  C 

06  C********************** *********************  **************************** 
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107 
‘  109 

;  io? 

5  <v 1 1  o 
'  m 

112 

113 

114 

115 

'  116 

117 

118 
'  119 

120 

121 

122 

123 

124 

125 

'  126 

127 

128 
12? 

130 

131 

132 
.  133 
.  :  134 

135 

136 

137 
13  3 
139 

1;  1 4  o 

141 

142 

143 

144 

145 

146 

147 

148 
14? 

150 

151 

152 

153 

154 

155 

156 

157 
153 

•-  159 


Now  get  the  inforciation  for  the  aperture  function  (  theta  , 
range,  and  lens  separation  ). 


Call  JCLR 
Call  Aptinf 


Multiply  every  other  element  of  the  two  functions  by  -1  to  force 
DC  terns  to  the  Middle. 


Call  JCLR 

Call  Invert  <  Ichose  ) 
C 

Cmmmm*  ****************** 

r 

C  Transforn  the  source. 

C 

C  *  * ******************* ****** 

c 

Call  FFTSrc  <  Ichose  ) 
Call  JCLR 


Multiply  the  transfer  ms  and  pupil  together. 


Cm****: 

c 


90 
9  0 
C 

c  *  *  * 
c 
c 
c 

Cm* 

C 


Write  <7 ,79) 

Fornatl/,  'Multiplying  pupil  by  source  F  F  T . ' ) 

Do  90  J  =  1  ,  256 

Do  80  I  =  1  ,  256 

RSourc  (  I  ,  J  )  =  RSourc  (  I  ,  J  )  *  RApert  (  I  ,  J  ) 
CSourc  (  I  ,  J  )  =  CSourc  <  I  ,  J  )  *  RApert  (  I  ,  J  ) 
Continue 
Continue 


***************************** 


Is  a  picture  of  this  product  needed? 


************************************* 


Write  (  7  ,  100  ) 

100  Fomat  (  /  ,$  ,  'Do  you  want  a  plot  of  the  product  of  the  FFT  of  the  so 
lurce  and  the  aperture',/, 'distributions  C  y / n  )  ?  ') 

Read  (  5  ,  60  )  Answer 
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Mon  Sep  23  09:22:15 

1985 

Wer  . 

4 . 02 

sy nap  t  .  f 

Page 

60 

If  ( 

Answer 

.  ne 

y  )  Goto 

105 

61 

Call 

PI tr s  t 

(  4 

,  0  ) 

62 

105 

Write 

(  7  , 

110 

) 

63 

110 

ForHat</;lj 

'Do 

you  wa  n  t  a 

hardcopy 

of  the  pr oduc  t  t  y/n 

1  ? 

64 

Read 

<  5  , 

60  ) 

Answer 

65 

If  ( 

Answer 

.  n  e 

y  )  Goto 

115 

66 

Call 

P 1 o  t  in 

67 

Call 

P 1 tr  s  t 

(  4 

,  0  ) 

68 

Call 

Plotof 

Inverse  F  F  T  this  resultant  Matrix  to  get  source  distribution 


115  Call  JCLR 

If  (  Ichose  .It.  5  )  Call  Invert  (  10  ) 
Call  IFTSrc 
Call  JCLR 


Write  (  7  ,  120  ) 

ForHat(/,$,'Do  you  desire  to  try  another  source  and  aperture  t  y/n 
13?') 

Read  (  5  ,  60  )  Answer 
If  (  Answer  .ne.  y  )  Goto  140 
Do  130  J  =  1  ,  256 

Real  (  J  )  =  0 
RiMag  (  J  )  =  0. 


130  I 

s 

1 

,  256 

R  S  o  u  r  c 

( 

I 

,  J  ) 

=  0 

CSo  urc 

( 

I 

,  J  ) 

=  0 

R  Ap  er  t 

( 

I 

,  J  > 

=  0 

Continue 
Lower  =  0 
Goto  10 
Call  Grphof 
Stop 
End 


0  Errors  detected 
!04  Source  lines  read 
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Subroutine  Grphin 


t  U  M  »  M  I  1 1  *  M  It  M  *  m  M  )  *  m  M  1 1 


This  connon  area  contains  the  required  data 


CoMMon  /  A  r  g  s  /  Real  <  256  )  ,  R  i  m  a  g  (  256  )  ,  RSourc  (  256  ,  25 
16  )  ,  CSourc  (  256  ,  256  >  ,  RApert  (  256  ,  256  )  , 

2  Plot  (  256  ,  256  )  ,  Ar  (  2  )  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 
C 

C 

C  Initialize  graphics  teminal  with  AGP  stuff  (  see  AGP  graphics 

C  Manual  ) . 

C 

c 

Call  JBEGN 

Call  JDINT  (1,3,  3Hwsp  ,  8  ,  8H/dev/tty  ,  0  ) 

Call  JWON  (  1  ) 

C 

C*« ********************************************************************* 

C 

C  Set  aspect  ratio. 

C 

C 

Call  JIMS  (  1  ,  254  ,0,0,2,  Idufi  ,  I  dun  ,  Ar  > 

Call  J  ASP  K  (  1.0  ,  Ar  (  2  )  ) 

C 

Ct****»t* ******  ********************************  ************************* 

c 

C  Set  up  viewing  references  and  rotate  axes. 

C 

C  *  ********************************  ************************************** 

c 

Call  JVDIS  (  1.0  ) 

Call  JPROJ  (  0  ,  0.3  ,  0.3  ,  -1.0  ) 

Wind  =  1.0 

Call  JWIND  <  -Wind  ,  Wind  ,  -Mind  *  Ar  (  2  )  ,  Wind  *  Ar  (  2  )  ) 


C ******************************* 

c 

C  Return  to  calling  prograti. 

C 

C******** *********************** 

C 

Return 

End 


********************************* 


****************************** 


ro  r.j  ro  ro 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 
19 

0 

1 

2 

3 


Subroutine  Twopnt 


C 

C  *  *  * 

c 
c 
c 
c 
c 

C  *  *  *  : 

c 


'***************************** 


This  subroutine  creates  a  two  point  source  by  Making 

RSourc  (112,128)  and  RSourc  (144,128)  both  1.0  and  leaving  the 

rest  of  the  array  0.0.  The  phase  is  also  0.0  at  all  points. 

I******************* 


CoMMon  /  Args  /  Real  (  256  )  ,  Rir-iag  (  256  )  ,  RSourc  (  256  ,  25 

16  )  ,  CSourc  (  256  ,  256  )  ,  RApert  (  256  ,  256  )  , 

2  Plot  (  256  ,  256  )  ,  Ar  (  2  )  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 

Write  (  7  ,  10  ) 

10  F or Ma t (/,*,' Enter  the  separation  distance  (  <  128  )  as  nnn  :  ') 
Read  (  5  ,  20  )  Length 
20  F  o  r  Ma  t  < 1 3 ) 

Ipos  =  128  -  INT  (  Length  /  2  ) 

RSourc  (  Ipos  ,  128  >  =  1.0 

RSourc  (  Ipos  +  Length  ,  128  )  =  1.0 

Return 

End 


0  Errors  detected 
23  Source  lines  read 


J»"*XTwtr*.v  urn*: m+.r^-+T*n.nm^K.  mu  mu  mu  m~ 
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Subroutine  Edge 


3  C  *  *  4 
A  C 

5  C 

6  C 

7  C 

8  C 

9  C  i  *  * 

10  c 

1 1 
12 

13 

14  C 

15  C  *  *  *  : 


This  subroutine  creates  an  edge  as  a  source.  This  is  done  by 
first  finding  how  wide  the  edge  is.  The  256  rows  X  Iwide 
coluMns  are  then  set  to  1. 


CoMMon  /  A  r  g  s  /  Real  <  256  >  ,  R  i  m  a  g  (  256  )  ,  RSourc  <  256  ,  25 

16  )  ,  CSourc  (  256  ,  256  )  ,  RApert  (  256  ,  256  )  , 

2  Flot  (  256  ,  256  )  ,  Ar  <  2  )  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 


17  C  Find  out  how  wide  the  edge  is  supposed  to  be. 

18  C 

19  C*i* *******************************************  ************************* 

20  C 

21  Write  <  7  ,  10  ) 

22  10  For  iia  t  <  /  ,  '  How  wide  is  the  edge  <  <  128  )  ?'  ,/,1,'Enter  answer  as  w 

23  lww  :  '  ) 

24  Read  <  5  ,  20  )  Iwide 

25  20  F orna  t ( 13  ) 

26  Do  40  J  =  1  ,  Iwide 

27  Do  30  I  =  1  ,  256 

23  RSo  urc  <  I  ,  J  )  =  1  .  0 

29  30  Continue 

30  40  Continue 

31  Length  =  256 

32  Return 

33  End 


0  Errors  detected 
33  Source  lines  read 
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Subroutine  Slit 


This  subroutine  creates  a  slit  for  a  source.  This  is  done  by 
first  obtaining  the  length  and  width  from  the  operator. 

The  Middle  Length  and  Iwide  rows  X  columns  are  then  set  to  i. 


Common  /  A  r  g  s  /  Real  (  256  )  ,  Rimag  (  256  )  ,  RSourc  (  256  ,  25 

16  )  ,  CSourc  (  256  ,  256  )  ,  RApert  (  256  ,  256  >  , 

2  Plot  (  256  ,  256  )  ,  Ar  (  2  >  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 
C 

C  t  n  I  M  H  M  1 1 1  «  I  M  M  m  M  1 1  t  M  M  M«  M  M  M  M  *  M  t  M  M  M  *  M  t  (  «  M  t  M  t  It  t  M  t  ( 

C 

C  Find  out  the  dimensions  of  the  slit. 

C 

C 

Urite  (  7  ,  10  ) 

10  For«at(/,  'What  is  the  length  and  width  of  the  slit  (  both  <  128  ) 

1  ,  Enter  as  111  www  :  '  ) 

Read  (  5  ,  20  )  Length  ,  Iwide 
20  Fomat(I3,lx,I3) 

C 

c 

C  Center  the  slit  around  the  point  128,128.  The  length  runs  along 

C  the  azimuth  and  the  width  along  the  range. 

C 

c 

J start  =  128  -  Iwide  /  2 
Istart  =  128  -  Length  /  2 
Do  40  J  =  Jstart  ,  128  +  J start 
Do  .3  0  I  =  Istart  ,  128  +  Istart 
RSourc  (  I  ,  J  )  =  1.0 
CSourc  (  I  ,  J  )  =  0. 

30  Continue 
40  Continue 
Return 
End 


0  Errors  detected 
43  Source  lines  read 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1  1 
12 

13 

14 

15 

16 
17 
19 

19 

20 
21 
22 
23 

2  4 

25 

26 

27 

28 

29 

30 

31 

32 

33 

3  4 

35 

36 

37 

38 
3? 

4  0 

41 

42 


C 

C 

C 

C 

C 

C 

C  *  *  t 
C 


c 

c 

c 

c 

c  *  * 
c 

10 


c 

C  * *  * 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C  *  *  * 

c 

c 


Subroutine  Circle 


This  subroutine  for ms  a  circular  source  in  the  RSourc  array. 

The  subroutine  utilizes  a  "shading"  routine  in  order  to  owercoMe 
the  rough  edges  caused  by  representing  circle  in  a  rectangular 
array  of  points.  See  cofiMents  later  for  details. 


CoMMon  /  Args  /  Real  (  256  )  ,  RiMag  (  256  )  ,  RSourc  (  256  »  25 

16  )  ,  CSourc  <  256  ,  256  )  ,  RApert  (  256  ,  256  )  , 

2  Plot  (  256  ,  256  >  ,  Ar  <  2  )  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 


Pronpt  operator  for  radius  Radius.  This  is  the  distance  out  frow 
the  center  of  the  arrays  that  the  circle  will  enconpass . 


Write  <  7  ,  10  ) 

For  fiat!/,*,'  Hhat  is  the  radius  of  the  source  (  <  65  )  t  ff.ff  1  ? 

1  '  ) 

Read  (  5  ,  20  )  Radius 
For«at<F6.2) 

Iupper  =  INT  (  Radius  ) 

******************************  ******  H  t  i  I  M  t  H  I  »  ♦  H  ♦  M  *  M  *********** 
The  interior  of  the  circle  is  now  filled  in. 

The  first  do  loops  deter  nine  how  far  out  froM  the  center  of  the 
array  (RSourc  (128,128))  should  be  filled  with  ones.  This  is 
deternined  to  be  the  point  just  prior  to  being  right  next  to 
the  outer  edge  of  the  circle.  The  conplex  part  is  assuned 
to  be  0.0. 

******************************************************************** 


4  3 

Do  70  I  =  0 

,  Iupper 

44 

Do  30  J  = 

0  ,  Iupper 

-  1 

45 

T  e  m  p  = 

I  *  I  +  J  * 

J 

46 

If  (  SQRT  (  Tenp  > 

.  g  t 

Radi 

47 

30 

Continue 

48 

40 

Do  50  K  = 

123  ,  128 

+  J 

4? 

RSourc 

(  128  +  I  , 

K  ) 

=  10 

5  0 

RSourc 

(  128  -  I  , 

K  ) 

=  1.0 

51 

50 

Continue 

52 

Do  60  K  = 

128  ,  128 

-  J  . 

,  -  1 

53 

RSourc 

(  128  +  I  , 

K  ) 

=  1.0 

For  tran/9000 
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54  RSourc  (  128  -  I  ,  K  )  =  1  .  0 

55  60  Continue 

56  70  Continue 

57  C***t *******************************************  ************************ 

58  C 

59  C  This  is  where  the  shading  begins.  The  the  distance  to  the  next 

60  C  column  is  still  short  of  the  radius  but  the  column  after  that  is 

61  C  greater  than  the  radius.  Shading  is  used  here  to  get  a  more 

62  C  accurate  representation  of  the  circle.  This  is  done  in  the 

63  C  following  Manner. 

64  C 

65  C  The  slope  of  a  line  running  from  the  last  point  filled  in  to  the 

66  C  first  point  outside  the  circle  such  that  the  line  crosses  the 

67  C  point  of  the  radius  of  the  circle  at  an  amplitude  of  0.5  is 

68  C  determined.  The  values  of  all  other  points  outside  the  radius 

69  C  is  extrapolated  froM  this  slope. 

70  C 


71 

72 


C 


73 

Do  100  I  =  0  ,  I  upper 

+ 

1 

74 

Value  =  1 

75 

Angle  =  ASIN  ( 

Float 

(  I 

) 

/  R  a  d  i  u 

3 

) 

76 

Trad  =  Radius  * 

COS 

( 

Angle 

) 

77 

Do  95  Ncols  =  1 

f 

I  u  p  p  e  r  + 

4 

79 

If  <  RSourc  (  128 

+ 

I 

,  129 

+  N  c  o  1 

3 

) 

gt . 

7  9 

Run  =  Abs  (  Float 

( 

Ncols  - 

1 

)  -  Trad 

> 

so 

Value  =  Value  -  ( 

0 

.  5 

* 

Run 

> 

81 

If  (  Value  .le. 

0 

.  ) 

Goto 

100 

92 

RSourc  (  128 

+ 

I 

> 

128 

4- 

Ncols 

) 

= 

Value 

93 

RSourc  (  128 

- 

I 

t 

128 

+ 

Ncols 

) 

s 

Value 

94 

RSourc  (  128 

+ 

I 

f 

128 

- 

Ncols 

> 

= 

Value 

85 

RSourc  (  128 

- 

I 

* 

129 

- 

Ncols 

> 

= 

Value 

96 

95  Continue 

87 

100  Continue 

89 

Re  turn 

89 

End 

Goto  95 


0  Errors  detected 
89  Source  lines  read 
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1 

2 

3 

4 

5 

6 

n 

/ 

9 

9 

10 
1  1 
12 

13 

14 

15 

16 
17 
13 
1? 
20 
21 
1  "3 

<_  i- 

23 

24 

25 

26 
27 
23 
2? 
30 
3  1 
3  2 
3  3 

34 

35 


C 


Subroutine  Other 


C 

C  This  subroutine  allows  the  operator  to  input  a  non-circular 

C  distribution  for  the  source. 

C 


CtMtMMMM 


c 

Co  mm  on  /  Args  /  Real  <  256  )  ,  R  i  m  a  g  (  256  )  ,  RSourc  (  256  ,  25 

16  )  ,  CSourc  (  256  ,  256  )  ,  RApert  (  256  ,  256  )  , 

2  Plot  (  256  ,  256  )  ,  Ar  <  2  )  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 
C 

C ********************************  *************************************** 

c 

C  The  operator  Must  input  the  values  of  both  the  real  and  i  m  a  g  i  n  - 

C  ary  parts  of  the  source  one  point  at  a  tine.  The  order  to  put 

C  the  points  in  is  one  row  at  a  tiMe. 

C 


C  *  *  *  *************************************************  ******************* 

c 

Write  <  7  ,  10  ) 

10  Forwatt 'Input  the  values  of  your  source.  You  have  a  256  X  256  are 
1 3 .  Enter  values  as  real  part,  in agin  ary  part  in  the  forwat  f.fff 
2 f f f .  Enter  values  by  row  starting  at  1,1  1,2  etc.') 

Do  40  I  =  1  ,  256 

Do  30  J  =  1  ,  256 

Read  <  5  ,  20  )  Rpart  ,  Cpart 
20  Format <F5 . 3  ,  lx  ,FS  .  3) 

RSourc  <  I  ,  J  )  =  Rpart 
CSourc  (  I  ,  J  >  =  Cpart 

30  Continue 

40  Continue 
Return 
End 


0  Errors  detected 
35  Source  lines  read 
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1  Subroutine  Plotin 


3 

4 

5 

6 

T 

/ 

8 

9 

10 
1 1 
12 

13 

14 

15 

16 
17 
13 
1? 
20 
21 
22 

23 

24 

25 


C  *****  * 

c 

C  This  coMMon  area  contains  the  required  data. 

C 

C*  **************************************  ******************************** 

C 

CoriMon  /  A  r  g  s  /  Real  <  256  )  ,  R i m a g  (  256  )  ,  RSourc  (  256  ,  25 
16  )  ,  CSourc  <  256  ,  256  )  ,  RApert  (  256  ,  256  )  , 

2  Plot  (  256  ,  256  )  ,  Ar  (  2  >  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 
C 

C**t  **************************************  ****************************** 

C 

C  Initialize  graphics  teminal  with  AGP  stuff. 

C 

c 

Call  Grphof 
Call  JBEGN 
Call  JDINT  (2,8 
Call  JWON  <  2  ) 

C 

c** ********************* 


8Hwsp . 7550  ,  13  ,  1 3H/d e v /p 1 t 755 0a  ,  0  ) 


************************************** 


26  C  Set  aspect  ratio. 

27  C 

28  C***t*  ***********************************  ******************************* 

2  9  C 

30  Call  JINS  (  2  ,  254  ,0,0,2,  Idu«  ,  Idurt  ,  Ar  ) 

31  Call  J  ASP  K  (  1 .  0  ,  Ar  <  2  )  ) 

32  C 

33  C**  ***********************************  ********************************** 

3  4  C 

35  C  Set  up  viewing  references  and  rotate  axes. 

36  C 

3  7  C  *  **********************************************  ************************ 

33  C 


3? 

Call 

JVDIS 

( 

1.0  > 

4  0 

Cal  1 

JPROJ 

( 

0  ,  0.3  ,  0.3  ,  -10  ) 

4  1 

W  l  n  d 

=  1.0 

42 

Call 

J  M  I  M  D 

( 

-Wind  ,  Wind  ,  -Wind  *  Ar  t  2  )  ,  Wind  *  Ar  (  2  )  ) 

43 

Call 

JNEWF 

44  C 

45  CM****  *********************************  ******************************** 

46  C 

47  C  Return  to  calling  p  r  o  g  r  a  m  . 

4  3  C 

49  Cl* *************************************** *♦♦********************'****** 

5  0  C 

51  Return 


•  I 
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5  C  This  subroutine  plots  the  contents  of  the  arrays  RSourc  and 

6  C  CSourc  on  the  HP  2623A  terminal.  If  the  aperture  function  is 

7  C  desired,  R Apert  is  plotted  instead.  These  arrays  contain  the 

9  C  data  for  either  the  source  distribution,  aperture  distribution, 

9  C  their  FFT  ,  the  product  of  their  FFTs,  or  the  final  obtained 

10  C  source  distribution.  An  appropriate  title  will  appear  for  each 

11  C  and  is  indicated  by  Iplot. 

12  C 

13  C  The  following  integer  arrays  are  used  to  pass  characters  to  the 

14  C  graphics  subroutine  JTEXM  which  writes  their  contents  to  the 

15  C  current  graphics  output  device.  AGP  requires  that  when  arrays 

16  C  are  to  be  written,  the  data  should  be  stored  as  INTEGER  *  2. 

17  C  The  data  in  the  data  statenents  is  therefore  in  Hollerith 

19  C  notation  for  this  reason.  The  progran  Manipulates  data  as 

19  C  nor  Mai  characters  through  the  use  of  character  arrays  which  are 

20  C  equivalensced  below  to  the  appropriate  INTEGER  *  2  array. 

21  C 

22  C  *  *  ******** * *************************************  *  *  *  ******************* * 

23  C 

24  Integer  *  2  Const  <  3  )  ,  Idelta  (  7  )  ,  Itheta  (  7  )  ,  Blank  , 

25  1  Marker  ,  Icount 

l:  to  w 

27  C*t***» ********************************************  ********************* 

23  C 

29  C  These  character  arrays  and  variables  are  used  to  plot 

30  C  character  strings  to  whatever  device  is  being  used  to  plot. 

31  C  CnorM  contains  the  (in  character  for m at)  nornalzation  constant 

32  C  R  n  o  r  m  .  Theta  contains  the  high  and  low  Units  on  the  theta 

33  C  input  by  the  operator  in  Aptinf.  InuM  contains  the  frequencies 

34  C  that  are  plotted  on  plots  of  the  aperture.  Dot  and  the  nuMbers 

35  C  (One,  Two,  etc.)  are  used  to  fill  the  above  arrays.  Answer  and 

36  C  y  are  used  in  deter  Mining  if  the  user  wants  another  plot. 

3  7  C 

38  C  *  ********************************  ************************************** 

39  C 

40  Character  CnorM  (  3  )  ,  Theta  (  14  )  ,  InuM  <  2  )  ,  Dot  ,  Two  , 

41  1  Three  ,  Four  ,  Six  ,  Eight  ,  One  ,  Answer  ,  y 

42  C 

43  C  ******  *  *****************************  *********************************** 

4  4  C 

45  C  These  are  the  Equivalence  state Ments  referenced  above.  This  a  1  - 

46  C  lows  the  data  to  be  Manipulated  in  the  progran  as  ASCII  but 

47  C  stored  in  INTEGER  *  2  arrays  for  use  by  the  JTEXM  subroutine. 

48  C 

4  9  C*M  *****************************************  *************************** 
3  0  C 

31  Equivalence  <  Const  (  5  )  ,  CnorM  (  1  )  )  ,  (  Const  (  6  )  , 

32  1  CnorM  (  3  )  )  ,(Const(7),Cnorn(5)),(Const<8), 

33  2  CnorM  <  7  )  )  ,  (  Itheta  (  1  )  ,  Theta  <!>),<  Itheta  <  2  )  , 
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54 

55 

56 

57 

58 

59 

60 
6  1 
62 

63 

64 

65 
6  6 
67 

69 
6? 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 
3  3 
84 
3  5 
86 
87 
39 
89 
9  0 
91 
9  2 


3  Theta  (  3  )  )  ,  <  Itheta  (  3  ) 

4  Theta  (7)  )  ,  <  Itheta  (5) 

5  (  Icount  ,  Inun  <  1  >  ) 


Theta  (  5  )  ) 
Theta  <  9  )  ) 


(  Itheta  (  4  ) 


25 


Co  mi  on  /  A  r  g  s  /  Real  (  256  )  ,  R  i  m  a  g  (  256  )  .  RSourc  (  256 

16  >  ,  CSourc  (  256  ,  256  )  ,  RApert  (  256  ,  256  )  , 

2  Plot  (  256  ,  256  )  ,  Ar  (  2  >  ,  Lower  ,  Iupper  ,  Length  ,  Iuide 

Data  B  1  a  n  k /2H  /  ,  Dot/'.'/  ,  Const  /2HR  n  ,2Hor  ,  2Hm  ,2H  =  , 

12H  ,2H  ,2H  ,2H  /  ,  I d e 1 1 a /2HDe  ,  2H 1 1 , 2Ha  /2HTh , 2He t , 2Ha  ,2H=  / 

1,  I theta/2H  ,2H  ,2Hto,2H  ,2H  ,2H  D,2Heg/  ,  I n u m/ '  1  '  ,  '  6 ' / , 

1  Two  /  '2 '/  ,  Three  / '3'/  ,  Four / ' 4  '/  ,  Six/ '6'/  ,  Eight/' 8'/  , 

1  One/'l'/  ,  y/'y'/ 


C 

C  1 1  * 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c  t  n  t  * 
c 


l 

c 

C  I 
c 
c 
c 

c  *  *  * 
c 


Determine  plot  title  as  follows; 

Iplot  =  1  Source  Irradiance  Distribution 
Iplot  =  2  Aperture  TransMi t  fanes  Function 
Iplot  =  3  FFT  of  Source  Distribution 

Iplot  =  4  FFT  of  Product  Source  and  Aperture  Distribution 
Iplot  =  5  Sanpled  Source  Distribution 

Initialize  Const. 


Do  1  I  =  5  ,  8 

Const  (  I  >  =  Blank 
Continue 

If  (  Iplot  .ne.  2  )  Goto  19 

!***t$***t**t«ttt*t«*.tt**tt***tttttttt*tt«*M***Mt*t*M*tt*t*ttt«* 
Put  Lower  and  Iupper  into  characer  for  Mat  for  plotting. 


9  3 

L  s  a  v  e 

Lower 

94 

Isa  ve 

= 

Iupper 

95 

Do  10 

I 

=  1,3 

96 

1 1  e np 

= 

Lower  / 

(  10 

*  *  (  3  -  I  ) 

) 

97 

I  c  h  a  r 

= 

I  teMp  + 

48 

98 

Theta 

( 

I  )  =  CHAR  ( 

Ichar  ) 

99 

Lower 

= 

Lower  - 

1 1  e  m  p 

*  (  10  *  *  ( 

3  - 

100 

If  ( 

Lower  .It. 

0  )  Lower  =  0 

10  1 

I  te Mp 

= 

Iupper  / 

(  10 

*  *  (  3  -  I  > 

) 

1  0  2 

I  c  h  a  r 

= 

1 1  e  Mp  + 

48 

103 

Theta 

( 

7  +  I  ) 

=  CHAR  (  Ichar  > 

104 

Iupper  = 

Iupper 

-  Itenp  ♦  <  10  ** 

(  3 

i  0  5 

If  ( 

Iupper  .It. 

0  ) 

Iupper  =  0 

"  06 

i  0 

Continue 

-  I  )  ) 


-  I  )  ) 
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Lower  =  Lsa ve 
I  upper  =  Isa ve 
Goto  85 


Find  Magnitude  of  input  data  and  find  the  normalization 
constant  Rnorn.  This  section  is  skipped  for  all  source  and 
aperture  plots  since  the  data  is  already  normalized  to  1. 


If  <  Iplot  .  le .  2  )  Goto  85 
Rnorn  =  0 . 

Do  30  J  =  128  ,  194 
Do  20  I  =  128  ,  194 

Plot  <  I  ,  J  )  =  <  RSourc  (  I  ,  J  )  *  RSourc  (  I 
1  CSourc  (  I  ,  J  )  »  CSourc  <  I  ,  J  )  >  *  »  0.5 

If  (  Plot  (  I  ,  J  )  .gt.  Rnorm  )  Rnorn  =  Plot  (  I 
Continue 
Continue 

Do  35  J  =  128  ,  194 
Do  35  I  =  128  ,  194 


J  )  + 


Plot  (  I 


J  )  =  Plot  (  I 


J  )  /  Rnorn 


The  following  lines  of  code  are  commented  out  for  now  as  per 
Maj  Mill's  instructions  .  This  code,  when  executed,  will  put 
the  data  in  dB  form  from  -100  dB  to  OdB.  All  data  less  than 
-100  dB  is  stored  as  -100  dB. 


If  <  Plot  (  I  ,  J  )  .ne.  0.  )  Then 

Plot  <  I  ,  J)  =20.  *  ALOG  (  Plot  <  I  ,  J  )  ) 

If  <  Plot  (  I  ,  J  )  .It.  -  100.  )  Plot  <  I  ,  J  )  = 
Else 

Plot  (I  ,J)=-100. 

Endif 

Continue 


Put  normalization  constant  in  character  form  for  writing  to 
screen  . 


Ikount  =  0 
Rsave  =  rnorm 
Temp  =  RnorM 
If  <  Temp  ,lt 


)  Goto  50 
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60 

61 

62 

63  50 

64 

65 

66 
67 


Ten p  =  T  e  m  p  /  10. 

Ikount  =  Ikount  +  1 
Goto  40 
lx  =  1 

Rnorn  =  RnorM  *  10000. 

Do  80  J  =  1  ,  Ikount  +  3 
If  (  Ikount  .gt.  0  )  Then 
If  (  J  .ne.  Ikount  )  Goto  60 


68 

69 

70 

71 

72 

73 

74 

75  60 

76 

77  70 

78 

79 

80 
81 
82 

83 

84 
35 
86 

87 

88 
99 

90 

91  80 

92  C 


Icon  =  46 
Goto  70 
Else 

If  (  J  .ne.  1  )  Goto  60 
I c  o  n  =  46 
Goto  70 
En  d  i  f 

Icon  =  RnorM  /  (  10.  <  *  <  Ikount  +  4  -  J  )  ) 

Icon  =  Icon  +  48 

If  (  Icon  .eq.  46  )  lx  =  lx  +  1 

Cnorn  (  lx  )  =  CHAR  <  Icon  ) 

If  (  Icon  .eq.  46  )  Then 
lx  =  lx  -  1 
Goto  60 
Else 

If  (  CnorM  (  lx  +  1  )  .eq.  Dot  )  Then 
lx  =  lx  +  2 
Else 

lx  =  lx  +  1 
En  d  i  f 
End  i  f 

If  <  Icon  . ne .  46  ) 

1  RnorM  =  Rnorn  -  (  Icon  -  48  )  *  (  10.  *  *  (  Ikount  +  4 

Continue 


J  )  ) 


93 

94 


95 

c 

Plot 

the  headings 

and  axes 

a  s 

appropriate . 

96 

c 

97 

c 

Ask 

the  user  if  he 

wants  2 

o  r 

three  diMensional  plots  and  tell 

98 

c 

AGP 

to  display  the 

data  as 

i  t 

is  called  for  here . 

99 


00 

c *************** 

0  1 

c 

02 

85 

Call  JIYQN 

0  3 

Write  (  7 

04 

82 

F  or  na  t  < /  ,  ' 

05 

13  ;  '  ) 

06 

Read  (  5  , 

0  7 

83 

ForMat(Il) 

03 

C 

0  9 

C  *  * 

************* 

10 

c 

1 1 


82  ) 


***************************************** 
Choose  a  window  size  for  the  plots, 
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213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 
2  25 
226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 
2  37 

238 

239 

240 

241 
2  42 

243 

244 

245 

246 

247 

248 

249 
25  0 

251 

252 

253 
25  4 

255 

256 


600 


610 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c  *  *  * 

c 

c 

c 

c 


Write  (  7  ,  600  ) 

Fornat(/,f,  'Choose  a  window  size  (default  is  1.4).  Enter  as  x.x 
1  '  ) 

Read  (  5  ,  610  )  Wind 
ForMat(F3.1) 

If  (  Wind  . le .  1 . E-7  )  Wind  =  1.4 

Call  JWIND  (  -Wind  ,  Wind  ,  -Wind  Mr  (  2  )  ,  Wind  *  Ar  (  2  >  > 


Set  the  pen  color.  This  is  ignored  for  the  2623a  since  it  is 
only  black  and  white.  The  plotter  has  eight  pens,  seven  of 
are  used  in  this  prograM.  They  a  are: 


1 

2 

3 

4 

5 

6 
*7 


Black 

□range 

Blue 

Light  Green 
Bark  Green 
Purple 
Red 


All  pens  are  3  mm  points  except  the  red  which  is  7  mm. 
JCOLR  is  the  graphics  call  that  changes  pen  color. 


Call  JCOLR  (  1  ) 
Call  J JUST  (  0.5 


0.0  ) 


If  (  Iplot  .  ne  .  2  )  Then 
Call  J 2 M O U  (  0.0 
Else 

Call  J2M0V  (  00 
End  i  f 


-0.8) 

-  0.62  *  Wind  /  1.4  ) 


The  above  statsnents  establish  the  pen  color  for  and  centers 
the  following  plot  headings  or  titles. 


257 

If 

( 

Iplot 

eq 

1 

) 

Call 

JTEXM 

< 

30 

258 

1  u  t  i  o  n 

) 

259 

If 

( 

Iplot 

eq  . 

2 

) 

Call 

JTEXM 

( 

31 

260 

1  n  c  t  i  o  n  ) 

261 

If 

< 

Iplot 

•  eq  ■ 

3 

) 

Cal  1 

JTEXM 

( 

26 

262 

In  ) 

263 

If 

( 

Iplot 

■  eq  • 

4 

) 

264 

ICa  1 1 

JTEXM 

<  33 

) 

33HFFT 

of  Source 

•96  5 

J 

If 

( 

Iplot 

eq  . 

5 

) 

Cal  1 

JTEXM 

( 

27 
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299  C 

290 

291  C 


295  C 
2  96  C  *  it 


303 

304 

305 

306 

307 

308  C 

309  C  »  t  ♦ 


ion  ) 


If  the  aperture  is  not  being  plotted/  print  the  normalization 
constant  and  label  the  axes  as  appropriate. 


If  (  Iplot  .ne.  2  )  Then 
Call  JCOLR  <  2  ) 

Call  J2M0V  (  0.5  *  Wind  /  1.4  ,  0.4  *  Wind  /  1.4  ) 
Call  JTEXM  (  9  ,  9HAmplitude  ) 

Call  J2M0V  (  0.5  *  Wind  /  1.4  ,  0.36  *  Wind  /  1.4  ) 
Call  JTEXM  (  10  ,  1 0 HNor m a  1 i z e d  ) 

Call  J2M0'J  (  0.4  «  Wind  /  1.4  ,  0.32  *  Wind  /  1.4  ) 


If  a  source  is  being  plotted,  do  not  print  Rnorn  since  the  data 
is  already  known  to  be  norwalized  by  default  anyways. 


If  (  Iplot  .ne.  1  )  Call  JTEXM  (  16  ,  Const  ) 


List  the  units  appropriately  for  V a  or  V  r  . 


If  (  <  Iplot  .It.  2  >  .or.  (  Iplot  .  g  t  .  4  )  )  Then 

Call  J2M0V  (  -  0.2  *  Wind  /  1.4  ,  0.45  »  Wind  /  1.4  ) 

Call  JCSIZ  <  0.015  ,  0.05  ,  1.0  ) 

Call  JTEXM  <  13  ,  13HVr  <Va>  is  in  ) 

Call  J2MQV  (  -  0.2  *  Wind  /  1.4  ,  0.41  ♦  Wind  /  1.4  ) 

Call  JTEXM  (  9  ,  9Hfiultiples  ) 

Call  J2M0V  (  -  0.2  *  Wind  /  1.4  ,  0.37  *  Wind  /  1.4  ) 
Call  JTEXM  (  13  ,  13Hof  500*pi/512  ) 

Call  JCSIZ  (  0  .  035  ,  0.05  ,  0.0  > 

Else 


tMtt**ttt**t***tttttttt«t*t** 


List  the  units  for  frequency. 


'.■■■  M 


« 
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319 

320 
\3  2  1 
'•322 

323 

324 

325 

326 

327 

328 


339 
34  0 

341 

342 

343 

344 
.345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 
37  0 
*71 


C 

c 

c 

c 


End  if 
Else 


When  plotting  the  aperture,  the  delta  theta  and  the  units 
must  be  listed. 


329 

Call 

JCOLR 

< 

1 

) 

330 

Call 

J3M0V 

( 

0  . 

35 

*  Wind 

I 

0  .  5 

331 

Call 

JCSI2 

( 

0  . 

01 

*  Wind 

/ 

1  .  2 

332 

Call 

JTEXM 

( 

14 

* 

Ide 1 ta 

) 

333 

Call 

J3M0V 

( 

0  . 

35 

*  Wind 

f 

0 . 45 

334 

Call 

JTEXM 

( 

14 

1 1  h  e  t  a 

> 

335 

Call 

JCSIZ 

( 

0  . 

015 

*  Wind 

/  1 . 2 

336 

Call 

J3MOV 

( 

- 

0  .  5 

*  Wind 

/  1  .  2 

337 

338 

Call 
End  if 

JTEXM 

( 

13 

t 

13HUn i t  s 

are 

Wind 
0  .  0.5 


1  .  2 
Wind 


♦  Wind  /  1.2 


0.0  ) 

/  1.2  , 

,0.0) 


0.05*  Wind  /  1 . 2 
0.5  *  Wind  /  1.2 


1.5  ) 


0.5  ) 

0.0  ) 


1 /m  ) 


C 

C 

C 

C 

c 

c 

c 


two 


Select  axis  projection  based  on  aperture  (  two 
di iiens.ons)  or  everything  else  (three  dimension).  When 
dimensional  plots  are  called  for,  the  three  dimensional 
projection  is  still  used  but  only  amplitude  and  azimuth  informa 
tion  are  plotted  along  the  azimuth  axis. 


700  If  (  Iplot  .  e  q  .  2 
1  Call  JPROJ  <  0 


0  .  0 


1  .  0 


-1.0) 


C******** 

c 

C  Draw  X 

C 

c 
c 

Ct*****t( 

c 


axis 


The  axis  is 


drawn  across  the  center  of  the  device 
for  the  aperture  plots  and  is  set  0.3  world  coordinate  units 
(see  AGP  graphics  manual  for  an  explanation  of  world  coordinate 


Call  JCOLR  (  3  ) 

If  (  Iplot  ,ne.  2  )  Then 
Call  J3H0V  (  -  1.0  ,  -  0,3 
Else 

Ca 1 1  J3M0V  (-10,0.0, 

End  if 

Call  JR  3  DR  (  2.0  ,  0.0  ,  0.0  ) 
Cal  1  J  JUST  (  0.0  ,  0.0  ) 

If  (  Iplot  ,ne.  2  )  Then 
Call  JCOLR  (  4  ) 


1.0  ) 


0.0  ) 
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*72 

(73 

574 

575 

576 
(77 
578 
5  7  9 
580 
181 
582 

183 

184 

185 
1 8  6 
>87 
>88 
>89 
>90 
191 

;  9  2 

193 

>94 

195 

196 

197 
;98 
;99 
00 

-401 
402 
4  03 

404 

405 
->0  6 

0  7 
08 
■09 
10 
11 
12 
-413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 


C 

C 

C 

C 

C 

C 

C  *  »  t 

c 


300 


Label  the  axes  for  anything  but  an  aperture.  If  plot  calls  for 
units  of  frequency,  label  theM  each  in  Multiples  of  32. 

If  units  of  V  are  called  for,  label  theM  in  eight  even 
increwents  froM  1  to  8. 


Do  84  I  =  1  , 
If  (  NUMplO 
Xpos  =  -  1.0 
Marker  =  10 
Call  J3MRK  ( 
Call  J3M0V  ( 
If  (  (  Iplot 
Goto  (  810 
Else 

Call  J3MRK 
If  (  NuMplo 
End  if 

Ypos  =  1.0  - 
If  (  NuMplo 
NuMplo  =  4 
Call  J3MRK  < 
Call  J3M0V  ( 
If  (  (  Iplot 
Goto  (  810 
Else 

Call  J3MRK 


8 

.  eq  .  4  )  Nump  1  o  =  3 
+  I  %  0.25 
+  I 

Xpos  ,-0.3  ,  1.0  ,2  ) 

Xpos  -  0.05  ,  -  0  35  ,  1.0  )  ' 

,gt.  2  )  .and.  <  Iplot  .It.  5  )  )  Then 
,  820  ,  830  ,  840  ,  850  ,  860  ,  870  ,  880  )  I 


Marker  ) 


(  Xpos  ,  -  0.35  ,  1.0 
.eq.  2  )  Goto  84 

I  »  0.25 

.eq.  4  )  Goto  84 


-1.0  ,  -  0  . 3  ,  Yp  os  ,  2  ) 

-  1.0  ,  -  0.35  ,  Ypos  +  0.05  ) 

.qt  .  2  )  .and.  (  Iplot  .It.  5  )  )  Then 

880  ) 


820 


830 


840 


(-1.0 


-0.35 


850 

Ypos 


860  ,  870 
Marker  ) 


Goto  84 
End  if 


310 

Call 

JTEXM  ( 

2  , 

2H3  2  ) 

If  ( 
Goto 

Nunplo 

800 

•  eq 

2  )  Goto 

84 

320 

Call 

JTEXM  < 

2  , 

2H64  ) 

If  < 
Goto 

NuMplo 

800 

eq  . 

2  )  Goto 

84 

330 

Call 

JTEXM  ( 

2  , 

2H96  ) 

If  ( 
Goto 

Nunplo 

800 

.  eq  . 

2  )  Goto 

84 

34  0 

Call 

JTEXM  ( 

3  , 

3H128  » 

If  ( 
Goto 

Nunplo 

800 

eq 

2  )  Goto 

84 

85  0 

Call 

JTEXM  < 

3  . 

3H160  ) 

If  < 
Goto 

Nunplo 

800 

eq  . 

2  )  Goto 

84 

860 

Call 

JTEXM  ( 

3  , 

3H192  ) 

I  f  < 
Goto 

NuMplo 

800 

eq 

2  )  Goto 

84 

870 

Call 

JTEXM  < 

3  , 

3H224  ) 

If  < 
Goto 

Nunplo 

800 

eq  . 

2  )  Goto 

84 
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125 

»26 

127 

»28 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 
16  0 
161 
162 

163 

164 

165 

166 
167 
163 

169 

170 

171 

172 

173 

174 

75 

76 


880 


84 


These  statenents  will  put  a  log  scale  on  the  plot  if  the  plot 
is  desired  in  dBs. 


Call  JTEXM  <3  ,  3H256  ) 

If  (  Nunplo  .eq.  2  )  Goto  84 
G  o  t  o  8  0  0 
Continue 

If  <  Nuwplo  .eq.  4  )  Nuwplo  =  3 
C 

C  *  *  * 

c 
c 
c 
c 

c  *  *  * 
c 
c 
c 
c 
c 
c 
c 

86 
C 

c  ********  t  *  *  t 
c 

Label  the  X  axis  fa  or  U  a  as  appropriate. 


If  (  Iplot  . It .  3  )  Goto  86 
Call  J3M0V  <  -  0.2  ,  0.77 
Call  JCOLR  (  6  ) 

Cal  1  JR3DR  (  0.4  ,  0.0  ,  0.0  ) 

Call  JCOLR  (  7  ) 

Call  JTEXM  (  5  ,  SH  3  dB  ) 

Call  J3MOV  (  0.98  ,  -  0.25 


0.0  ) 


1.0  > 


c 

c 

c  *  * 
c 


Else 


Cal  1 

J3M00 

(  0  . 

95  ,  -  0  .  1  , 

0.0  ) 

End  if 

[f  (  ( 

Iplot 

.  g  t  . 

2  )  .and.  ( 

Iplot 

Call 

JCOLR 

(  2 

) 

Call 

J3M0V 

(  0  . 

99  ,  -  0.25  , 

1.0  ) 

Call 

JCOLR 

(  1 

) 

Call 

JTEXM 

(  i 

,  1 H  f  ) 

Call 

J  JUST 

(  0  . 

0,0.5) 

Call 

JTEXM 

(  i 

,  1  Ha  ) 

•  lse 

Call 

JCOLR 

(  i 

) 

If  ( 

Iplot 

.  n  e  . 

2  )  Then 

Call 

JCSIZ 

(  0 

.060  ,  0.08  , 

0  .  05 

Call 

JTEXM 

(  2 

,  2HVa  ) 

Else 

Call 

JCSIZ 

<  0 

. 070  ,  0.10  , 

0.4  ) 

Call 

JTEXM 

(  2 

,  2Hf  a  ) 

Cal  1 

JCSIZ 

(  0 

.035  ,  0.05  , 

0.0  ) 

It.  5  )  )  Then 


End  if 
End  if 


Ct*M*t** ******************************  ********************************* 

C 

C  Draw  anplitude  axis  if  not  plotting  the  aperture. 


Fortran/9000 
Ve r  .  4.02 


Hon  Sep  23  09.‘35i3i  1985 
pltrst.f  Page  10 


473 

479 

490 

431 

432 
483 
434 
495 
436 
497 

433 
439 

490 

491 

492 

493 

494 

495 

496 

497 
493 
499 
50  0 

501 

502 

503 

504 

505 

506 

507 


527 

528 

529 

530 


C*  *  * 

c 

If  <  Iplot  ,ne.  2  )  Then 
Call  J3MOV  (-1.0  ,-0.3,10) 

Else 

Call  J3MOV  (  0.0  ,  0.0  ,  0.0  > 

En  d  i  f 

Call  JCOLR  (  3  ) 

Call  JR3DR  (  0.0  ,  1.0  ,  0.0  > 

C 

Cl  M  M  t  M  n 
C 

C  Label  the  anplitude  as  being  the  Modulus 

C 

C  ********** 

c 

If  (  Iplot  .ne.  2  )  Then 
Call  JCSIZ  (  0.06  ,  0.08  ,  0.05  ) 

Call  J3M0V  (  -  0.95  ,  0.72  ,  1.0  ) 

Call  JTEXM  <  7  ,  7HModulus  ) 

Else 

Continue 
En  d  i  f 


C 

C  *  *  * 
c 
c 
c 


Draw  the  Y  axis  only  for  aperture  and  three  diMensional  plots. 


508 

If  (  Iplot  . 

ne 

.  2 

) 

Then 

509 

Call  J3M0V 

( 

-  1 

.  0 

*  " 

0  .  3 

,  1.0 

510 

Else 

511 

Call  J3M0V 

( 

0  .  0 

l 

0.0  , 

-  1 

0  ) 

512 

Endif 

513 

Call  JCOLR 

< 

3  ) 

514 

If  (  Iplot  . 

eq 

2 

) 

Then 

515 

Call  JR3DR 

< 

0  .  0 

1 

0.0  , 

2 . 0 

) 

516 

Else 

517 

If  <  N  U  M p 1 0 

ne  . 

2 

)  Call  JR3DR  ( 

5  18 

Endif 

519 

If  (  Iplot  . 

ne 

.  2 

) 

Then 

Call  J  JUST 

( 

0  .  0 

* 

0.0  ) 

521 

Call  J3I10V 

( 

-  0  . 

35 

t 

0  35 

*  ” 

522 

Else 

523 

Call  JJUST 

( 

0  .  0 

I 

0.0  ) 

524 

Call  J3M0V 

( 

0.01 

* 

-  0  . 

in 

<=• 

1.0  ) 

525 

Endif 

526 

Call  JCOLR 

( 

1  ) 

0  .  0 


.  0  ) 


C ********  it** 

c 

C  Label  the  Y  axis  fr  or  Vr  as  appropriate. 


*********** **»**»»***»*»***»»*** 
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C 

c 

C  *  *  * 

c 


If  (  (  Iplot  .gt.  2  )  .and.  (  Iplot  .It.  5  >  )  Then 

If  (  N  u  m  p  1  o  .ne.  2  )  Then 
Call  JTEXM  (  1  ,  lHf  ) 

Call  J3MOV  <  -0.85  ,  -  0.36  ,  -  1.0  ) 

Call  JTEXM  (  2  ,  2H  r  ) 

Else 

Continue 
End  i  f 
Else 

If  <  Iplot  . eq.  2)  Then 
Call  JCSIZ  (  0 . 070  ,  0.10  ,  0.4  ) 

Call  JTEXM  (  3  ,  3H  fr  ) 

Call  JCSIZ  (  0 . 035  ,  0.05  ,  0.0  ) 

Else 

Call  JCSIZ  (  0.060  ,  0.08  ,  0.05  > 


0.4  ) 


0.0  ) 


0.05  > 


If  (  Nutiplo  .ne  ,  2  )  Call  JTEXM  (  3 
Call  JCSIZ  (  0  .  035  ,  0  .05  ,  0.0  ) 


3H  Ur  ) 


End  i  f 
E  n  d  i  f 


Plot  the  data. 


Xincr  =  1.  /  128. 

Yincr  =  1.  /  128. 

If  (  Iplot  .gt.  2  )  Goto  ISO 

Goto  (  90  ,  100  ,  110  ,  140  ,  150  ,  150  )  Ichose 


C 

C  t  H  M 


Plot  a  point  source 


********************** 


C 

CM* 


Call  JCOLR  <  1  ) 

Call  J3MRK  (-1.0, 0.7, 1.0, 3) 
Goto  500 


t $t t ******** t* * t t ****************** 


L 

CM* 


Plot  a  two  point  source. 


t **********************  ******************* 


Call  JCOLR  (  1  ) 

Ypos  =  10  -  (  Float  <  Length  >  /  2.  )  *  Vincr 
Call  J3MRK  (  -  1.0  ,  0  7  ,  Yd  os  ,  3  ) 


*v\ 
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Goto  500 
C 

C 

C  Plot  an  edge. 

C 

C 

110  Call  JCOLR  <  5  ) 

Do  120  J  =  1  ,  I  NT  (  Iwide  /  2  ) 
Xpos  =  -  1.0  +  J  *  Xincr 
Call  J3M0V  (  X  p  o  5  ,  0.7  ,  1.0  ) 

If  (  Nuoplo  .eq.  2  )  Goto  120 
Ca  1 1  J3DRW  <  Xpos  ,  0.7  ,  -  1.0  > 
120  Continue 

Call  JCOLR  (  1  ) 

Do  130  1=1,  256 
Ypos  =  1.0  -  I  *  Yincr 
Call  J3M0V  <  -  1.0  ,  0.7  ,  Ypos  ) 
Call  J3DRW  <  Xpos  ,  0.7  ,  Ypos  ) 

If  <  Nunplo  eq.  2  )  Goto  500 
130  Continue 
Goto  500 


C 

C  Flot  a  slit. 

C 

C 

140  Ypos  =  1.0 

Xpos  =  -  1.0 

Call  J3M0V  (  -  1.0  ,  0.7  ,  1.0  ) 

Dist  =  Float  (  Iwide  )  *  Xincr  /  2. 
Call  JCOLR  (  1  ) 

Do  145  I  =  1  ,  INT  (  Length  /  2  )  +  1 
Call  JR3DR  (  Dist  ,  0.0  ,  0.0  ) 

Ypos  =  1.0  -  I  *  Yincr 

If  <  NuMplo  .eq.  2  )  Goto  500 

Call  J  3  M  O  V  (  Xpos  ,  0.7  ,  Ypos  ) 

145  Continue 
Goto  500 


*'**,*******♦»«**»*»**»* 


Plot  a  circle 


C  ********  ************************** 

C 

150  Xincr  =  1.0  /  16. 

Yincr  =  1.0  /  16. 

Call  JCOLR  (  1  ) 

Do  170  J  =  112  ,  192 


84 
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37 

39 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
6  1 
62 

63 

64 

65 
6  6 
67 
69 

69 

70 

71 

7  7 

.  t_ 

73 

74 

75 

76 

77 

78 

79 
9  0 
81 
82 

83 

84 

95 

96 

97 
99 
89 


Do  160  I  =  112  , 

If  (  Iplot  .  eq  . 

If  <  R Apert  (  I  , 

Xpos  =  -  1.0  +  (  J  - 
Ypos  =  1.0  -  >.  I  - 

Call  J3HRK  (  Xpos 

If  (  I  . g  t .  147  )  Goto  170 

If  (  J  .gt.  147  )  Goto  160 

Else 

If  (  <  I  It.  128  )  .  or  .  ( 

If  (  RSourc  (  I  ,  J  )  eq. 

If  (  (  Iplot  .eq.  1  )  .and.  ( 
1)  )  Goto  170 

Yincr  *  1  /  32. 

Xincr  =1/32. 

Ypos  =  1 . 0  -  <  I  -  128  >  * 

Xpos  =  -  1 . 0  +  (  J  -  128  ) 

Zpos  =  RSourc  (  I  ,  J  >  - 
Call  JCSIZ  <  0.070  ,  0.10 
Call  J3MRK  (  Xpos  ,  Zpos  , 
End  if 

160  Continue 
170  Continue 

If  <  Iplot  ,ne.  2  )  Goto  500 
Call  JCOLR  (  4  ) 

Do  175  1=1,4 
Pos  =  0.25  *  I 
If  (  I  .eq.  1  )  Then 

InuM  (  1  )  =  One 
Inuri  (  2  )  =  Six 
Else 

If  (  I  .eq.  2  )  Then 
InuM  (  1  )  =  Three 
InuM  (  2  )  =  Two 
Else 

If  (  I  eq.  3  )  Then 
Inufi  (  1  )  =  Four 
In  urt  (  2  >  =  Eight 
Else 

InuM  <  1  )  =  Six 

In  um  (  2  )  =  Four 
En  d  i  f 
E  n  d  i  f 
E  n  d  i  f 
C 


0.0,1) 


J  It.  128  )  )  Goto  160 
0  )  Goto  160 

Nuttplo  eq.  2  )  and.  (  I  .gt.  128 


Yincr 
*  Xincr 
0  .  3 

,0.4) 

Ypos  ,  1  ) 


192 

2  )  Then 

J  )  eq.  0  )  Goto  160 
112  )  *  Xincr 
112  )  *  Yincr 
,  Ypos 


C  M  I  1  t  «  M  K  M  *  Mt  I  »  t  *  M  M  MH  Mil*  t  MM  »  t  M  t  *  M  M  t  H  M  «  H  M  M  M  I  M  M  t 

c 

C  Label  the  axes  for  an  aperture  in  Multiples  of  16. 

C 

C  M  M  t  MM  M  M  M  M  M  t  M  M  MM  M  «  M  *  M  MMMM  M  M  M  M  M  I  M  M  M  MM  M  M  M 

c 

Call  JCSIZ  (  0.060  ,  008  ,  0.01  > 

Call  J3MRK  (  Fos  ,  0.0  ,0.0,2) 
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Call  J3MRK  <  0.0  ,  0.0  ,  Pos  ,  2  ) 

Call  J3MRK  <  -  Pos  ,  0.0  ,  0.0  ,  2  ) 

Call  J3MRK  (0.0,  0.0, -Pos, 2) 

Call  J JUST  (  0.5  ,  0.0  ) 

Call  J3MOV  <  Pos  ,  0.03  ,  0.0  ) 

Call  JTEXM  <  2  ,  Icount  ) 

Call  J  JUST  (  0.0  ,  0.0  ) 

Call  J3MDV  <  -  0.15  ,  Pos  ,  00  ) 

Call  JTEXM  (  2  ,  Icount  ) 

175  Continue 

Call  JPROJ  (0,0. 3,  0.3,  -1.0) 

Goto  500 


C 


ISO  Xincr  =  1.0  /  32. 

Yincr  =  1.0  /  32. 

C  Zinc  =1/100. 

Call  JCOLR  <  1  ) 

Do  200  I  =  129  ,  193  ,  2 

Ypos  =  1.0  -  (  I  -  129  )  t  Yincr 

Z  =  <  Plot  <  I  ,  129  )  /  Rsave  )  -  0.3 

C  Z  =  1  +  Plot  <  I  ,  97  )  *  Zinc  -  0.3 

Call  J3MOU  (  -  1.0  ,  Z  ,  Ypos  > 

Do  190  J  =  129  ,  193  ,  2 
Xp  os  =  -  1.0  +  <  J  -  129  )  *  Xincr 

Z  =  (  Plot  (  I  ,  J  )  /  Rsave  >-0.3 

C  Z  =  1  +  Plot  (  I  ,  J  )  »  Zinc  -  0.3 

Call  J3DRW  (  Xpos  ,  Z  ,  Ypos  ) 

190  Continue 

If  (  Nunplo  .eq.  2  >  Goto  500 
200  Continue 
C 

Call  JCOLR  (  5  ) 

Do  220  J  =  129  ,  193  ,  2 
Xpos  =  -  1.0  +  (  J  -  129  )  t  Xincr 
Z  =  (  Plot  (  129  ,  J  >  /Rsave  )  -  0.3 
C  Z  =  1  +  Plot  (  129  ,  J  )  *  Zinc  -  0.3 

Call  J 3 M 0 V  <  Xpos  ,  Z  ,  1.0  > 

Do  210  I  =  129  ,  193  ,  2 
Ypos  =  1.0  -  (  I  -  129  )  *  Yincr 
Z  =  <  Plot  <  I  ,  J  )  /  Rsave  >-0.3 
C  Z  =  1  +  Plot  (  I  ,  J  >  *  Zinc  -  0.3 

Call  J3DRW  <  Xpos  ,  Z  ,  Ypos  > 

210  Continue 
220  Continue 
500  Call  JIVOF 

Call  JCOLR  (  1  ) 
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43 

Call  JWIND 

(  -Wind  ,  Wind 

,  -Wind  *  Ar  (  2  >  , 

Wind  » 

Ar  <  2  ) 

44 

Write  (  7  , 

510  ) 

45 

510 

Foma  t  <  /  ,  *  , 

'Would  you  like 

to  replot  t  y/n  J  ? 

'  ) 

46 

Read  (  5  , 

520  )  Answer 

47 

520 

Foma  t  <A1  ) 

48 

If  (  Answer 

.  eq  .  y  >  Goto 

85 

49 

Return 

50 

End 

0  Errors  detected 
750  Source  lines  read 
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1 

O 

i— 

3 

4 

5 

6 

9 


Subroutine  Plotof 


C 


C  *  t  * 


C 

C  This  subroutine  turns  off  the  graphics  when  the  Main  prograM 

C  is  f in  ished  . 

C 


c  *  *  *  *  «  * 


?  c 


i  0 

CofiMon  /  A  r  g  s 

/  Real 

<  256 

) 

,  RiMag  (  256  )  , 

RSo  urc  ( 

256  ,  25 

1 1 

16  )  ,  CSo  urc  < 

256  , 

256  ) 

» 

RApert  (  256  ,  256 

)  , 

12 

13  C 

2  Plot  (  256  , 

256  ) 

Ar  ( 

2 

>  ,  Lower  ,  I  upper 

,  Length 

,  Iwide 

14 

Call  JWOFF  (  2 

) 

IS 

Call  JWEND  <  2 

) 

16 

Call  JEND 

17 

Call  Grphin 

18 

Return 

19 

End 

0  Errors  detected 
19  Source  lines  read 
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Subroutine  Aptinf 


This  subroutine  foms  the  aperture  according  to  the  operator's 
desires.  The  operator  gives  the  range  of  theta  Measured. 

The  coMputer  then  determines  what  frequencies  are  sampled  to 
ensure  that  the  samples  fall  exactly  on  a  point  in  the  array 
and  not  somewhere  in  between.  This  is  done  because  the  FFT 
subroutine  cannot  work  in  radial  coordinates.  Only  frequencies 
that  lie  at  specific  points  in  the  aperture  are  measured  due 
to  the  bandpass  limited  nature  of  the  system. 

The  subroutine  determines  the  MaxiMUM  theta  based  on  the  lens 
separation  (  Sep  )  and  range  <  Range  ).  These  figures  are  input 
by  the  user.  The  theta  that  the  user  wishes  to  sample  within 
is  then  asked  for.  The  two  thetas  are  compared  to  see  that 


Vv\l 


that 

the  theta  input  by  the  user  falls  within  the  limits  of  the 
system  as  determined  by  the  range  and  separation.  The  spatial 
frequencies  to  be  sampled  are  based  on  filter  frequencies  that 
lie  between  8  and  12  umicrons.  The  highest  spatial  frequency 
is  64  1/m.  This  limit  was  arrived  at  by  considering  the 
best  Fourier  transform  vs  the  width  of  the  pupil.  A  pupil  width 
of  32  was  found  to  be  optimum.  64  1/m  was  the  highest  frequency 
that  could  be  sampled  due  to  the  bandpass  of  the  system  and 
assuming  a  range  of  no  less  than  1  Km.  This  frequency  cor¬ 
responds  to  the  outside  edge  of  the  pupil.  Therefore,  the  fre¬ 
quencies  measured  are  multiples  of  4  (  4  x  16  (  radius  )  =  64  ) . 

The  subroutine  determines  the  minimum  and  maximum  spatial 
frequencies  based  on  the  relationship 

f  =  (2  *  fm  /  c  )  *  sin  <  alpha  /  2  )  unit  vector  f 

where  fm  /  c  is  1  /  wavelength  and  alpha  is  twice  the  angle 
defined  by  the  range  from  a  point  half  way  between  the  lenses 
to  the  target  and  the  distance  from  the  target  to  a  lens. 

(See  ERIM  report.)  The  nearest  multiple  of  4  is  found  and  all 
points  within  the  specified  theta  that  lie  between  the  minimum 
and  maximum  spatial  frequency  are  set  to  one  indicating  that 
the  corresponding  frequency  is  sampled. 


********** 


********************* 


(  256 
256  ) 
Ar  ( 


Common  /  Args  /  Real 
16  )  ,  CSourc  <  256  , 

2  Plot  (  256  ,  256  ) 
Character  Answer  ,  y 
Data  y /  '  y  '  / 

Pi  =  3  1415926535 


tttt ***************************  ***** 


.  R  l  m  a g  (  256  )  ,  R  S  o  u  r  c  (  ; 
RApert  (  256  ,  256  )  , 

)  ,  Lower  ,  Iupper  ,  Length 


Iwide 


89 
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Get  the  range  and  lens  separation  distance  in  kiloneters  and 
neters  respectively  f  r  o  m  the  user. 


Lower  =  -  1 
Write  C  7  ,  2  ) 

ForMa t ( / ,  $  , 'En ter  range  (Km)  and  lens  separation  <m>  as  x.x  x.xx 
1  '  ) 

Read  (5,3)  Range  ,  Sep 
Forwa t (F3 .  1 ,  lx  ,F4 .2) 


Calculate  the  angular  separation  (  in  radians  )  of  the  lenses 
for  the  input  range  and  separation. 


A 1 pha2  =  ATAN  (  (  Sep  /  20  0  0.  >  /  Range  ) 


Calculate  the  m  a  x  i  m  u  m  theta  that  can  be  covered  assuMing  a 
collector  speed  and  stability  and  range  specified  by  the  user 
(  Dtheta  in  degrees  ) 


C 

c  *  *  * 

c 

c 

c 

CM* 

c 


Write  (  7  ,  4  ) 

For Mat (/,$,' Enter  collector  speed  (ft/s)  and  stability  as  xxxx.x  x 

1  x  .  x  !  '  ) 

Read  (  5  ,  5  )  Speed, Stabil 
ForMat(F6.1,lx,F41) 

Dtheta  =  360.  *  ATAN  (  (  Stabil  *  Speed  *  .0003048  /  2.  ) 

1  /  Range  )  /  Pi 


Get  the  user's  lower  and  upper  Units  on  theta. 


Write  (  7  ,  9  ) 

For  m  a  t ( /  ,  'Enter  theta  for  your  particular  aperture.') 

Write  (  7  ,  20  ) 

For*iat(/,  'What  range  does  theta  lie  within  froM  0  to  180  degrees?' 
1  >/,/,*,  'Enter  your  answer  as  nnn  nr.n  :  ') 

Read  (  5  ,  30  )  Lval  ,  Iupper 
ForMat(I3,lx,I3) 

If  (  (  Lval  It.  Lower  )  .or.  (  Lower  .  e  q  .  -  1  )  )  Lower  =  Lval 
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107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 
1  32 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 
143 

149 

150 

151 

152 

153 

154 

155 

156 

157 
159 
159 


Check  to  see  if  the  theta  input  by  the  user  falls  within  the 
calculated  theta  froti  above. 


34 


If  <  (  Iupper  -  Lval  )  .le.  Dtheta  )  Goto  35 

Write  (  7  ,  34  )  Dtheta 

Forrtat</,'You  have  exceeded  the  m  a  x  i  m  u  m  Dtheta  of  '  ,  F  7 . 3  ,  '  degrees 
1  for  the  range  e n t e r e d Tr y  again.') 

Goto  i 


C 
C 
C 
C 

C  *  *  t 

c 

35 


Check  to  see  that  theta  lies  between  0  and  180  degrees.  If  not, 
ask  the  operator  to  input  theta  again. 


40 


If  (  <  Lval  .It.  0  )  .or.  (  Lval  .  g  e  .  180  )  or.  (  Iupper  .le 
1  0  )  .or.  (  Lval  .ge.  Iupper  )  >  Then 
Write  (  7  ,  40  ) 

F or «a t (  '  L i« i ts  of  theta  are  out  of  bounds.  Try  again.') 
Goto  10 
Else 

Continue 
End  l  f 


C 


C  *  *  *  *  *  * 

C 

C  Convert  the  bounds  on  theta  to  radian  quantities 

C 

c  m  *  t  H  t  M  M  m  »  M  M  t 

c 

Radlow  =  Lval  *  3.1415927  /  160. 

Padhi  =  Iupper  *  3.1415927  /  ISO. 

C 

C **********************************************  ************************* 

C 

C  Find  the  mxiHUH  and  nininuH  spatial  frequencies  U p r a d  and 

C  rlorad .  All  points  equal  to  or  lying  between  these  two  points 

C  are  set  to  one.  rlwav  and  Uwav  are  the  lower  and  upper 

C  wavelengths  satipled  by  the  systen  . 

C 

Ctt** ************************************************ ******************* 

c 

rlwav  =  8E-6 
Uwav  =  1 2  E - 6 

Uprad  =  <  2.  /  rlwav  )  *  Sin  (  Alpha2  > 
rlorad  =  (  2.  /  Uwav  )  *  Sin  (  Alpha2  ) 

C 

c**(  **********************************  ************* 

c 
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Check  to  see  that  the  user  has  not  asked  for 
the  radius  of  the  pupil. 


frequency  outside 


ClUMMMttUK 

C 

If  (  Uprad 
Write  (  7 


If  (  Uprad  . le .  64  )  Goto  60 
Write  (  7  ,  50  ) 

Format! 'Your  combination  of  lens  separation  and  range  caused',/, 
i'the  upper  spatial  frequency  to  fall  outside  the  aperture.',/, 
2'Please  try  again.') 

Goto  1 


c  H  t  »  *  It  t  M  M  M  »  U  m  (  l  *  )  U  t  m  M  t  M  U  *  M  M  )  t  M  t  M  t  t  Mt  m  M  M  M  «  »  t  U 

c 

C  Calculate  the  inner  and  outer  radii  <  squared  )  between  which 

C  the  spatial  frequencies  to  be  sampled  lie.  The  maximum  radius 

C  is  rounded  up  and  the  lower  radius  rounded  down  to  give  the 

C  largest  n  u  fiber  of  frequencies  sanpled  possible.  This  gives  the 

C  benefit  of  the  doubt  in  the  system's  favor. 

C 

CimiiimmtMmttimtiMtimimmmmMtiimtMtmmmi 

c 

60  Const  =  ANINT  <  Uprad  /  4.  ) 

Const  =  Const  *  Const 
Other  =  AINT  (  rlorad  /  4.  ) 

Other  =  Other  *  Other 


C 

C  Mow  determine  where  the  lower  and  upper  bounds  of  theta  lie. 

C  The  possibilities  are  0  to  45  ,  45  to  135,  and  135  to  180 

C  degrees. 

C 

c 

C** *******************************************  ************************** 

c 

C  0  to  45  degrees.  First  compare  Radio w  (  the  lower  bound  on 

C  theta  )  to  the  angles  formed  by  the  inverse  tangent  of  the 

C  column  divided  by  the  radius  This  determines  whether  or  not 

C  the  section  for  0  to  45  degrees  should  be  carried  out.  If 

C  Radlou  fall  between  these  limits,  the  program  jumps  to  100  to 

C  determine  what  points  in  the  aperture  are  to  be  set  The  exact 

C  starting  point  is  determined  in  the  section  of  the  program 

C  starting  at  100. 

C 

C* ******************************************  **************************** 

c 

Do  70  J  =  0  ,  16 

Temp  =  SORT  (  J  *  J  /  Const  ) 

Angle  =  ATAN  (  Temp  ) 

If  (  Angle  .gt.  Radlow  )  Goto  100 
70  Continue 
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213 

214 
>.215 

2  1  6 

217 

218 

219 

220 
221 
2  22 

223 

224 

225 

226 
227 
229 

229 

230 

231 

232 

233 
23  4 

235 

236 

237 

238 
,%239 
V  240 

241 

242 

243 

244 

245 
2  46 

247 

248 

249 

250 
2  51 

252 

253 

25  4 

255 

256 

257 
259 
2  59 

26  0 
261 
2*2 
2*3 
2*4 
"65 


C 

C  45  to  135  degrees.  See  if  lower  bound  on  theta  falls  in  here. 

C  If  Radlow  was  greater  man  45  degrees,  the  subroutine  continued 

C  right  on  into  this  do  loop.  The  angle  is  conputed  as  above 

C  except  that  first  pi  Minus  the  angle  is  checked  and  then  pi  plus 

C  the  angle  is  checked.  The  subroutine  then  j  u  n  p  s  to  the  section 

C  where  these  angles  are  dealt  with  (  9  150  ). 

C 


c 


Do  80  12  =  16  ,  -  16  ,  -  i 
TeMp  =  SQRT  (  12  *  12  /  Const  ) 

If  (  12  .ge.  0  )  Then 
Angle  =  (  Pi  /  2.  )  -  ATAN  <  TeMp  > 

Else 

Angle  =  (  Pi  /  2.  )  +  ATAN  (  TeMp  ) 

Endif 

If  (  Angle  .gt.  Radlow  )  Goto  150 
SO  Continue 


C 

C* **********************************  ************************************ 

c 

C  While  the  lower  bound  Must  lie  in  here  by  process  of  elinina- 

C  tion  (  135  to  180  degrees  ),  this  check  Makes  sure  the  user 

C  did  not  goof  when  inputting  theta.  The  angle  is  found  as 

C  before  and  lies  between  135  and  180  degrees. 

C 

C*t*M ****************************************  ************************** 

c 


Do  90  J2  =  16  ,  0  ,  -  1 

TeMp  =  SQRT  (  J2  *  J2  /  Const  ) 

Angle  =  Pi  -  ATAN  (  Tewp  ) 

If  (  Angle  .gt.  Radlow  )  Goto  190 
90  Continue 
Goto  10 
C 

C»*«**** ************************************  ************************** 

C 

C  This  is  where  the  frequencies  that  will  be  allowed  to  pass  are 

C  deternined.  All  the  possible  slopes  t'roM  0  to  infinity  are 

C  calculated  and  any  point  lying  along  that  slope  that  is  not 

C  outside  the  outer  radius  or  less  than  the  inner  radius  of  the 

C  pupil  is  set  to  one  This  section  pertains  to  slopes  starting 

C  at  0  to  45  degs.  The  angle  of  the  slope  is  found  by  calculating 

C  the  inverse  tangent  of  the  coluMn  /  row 

C 

C  For  e  x  a  m  p  1  e  ,  the  first  check  is  to  be  sure  that  the  slope's 

C  angle  is  not  less  than  the  lower  or  greater  than  the  upper 

C  bound  on  theta.  This  deter  nines  at  what  angle  the  aperture 

C  will  begin.  If  this  condition  is  Met,  fill  in  all  points 

C  along  the  slope  greater  than  the  square  root  of  Other  and  less 
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2  66 
267 
263 
269 
27  0 
271 
^  72 

273 

274 

275 

276 


t  7  g 
0  7  g 
230 
23  1 

232 

233 

234 

235 

236 


238 

239 
29  0 

291 

292 

293 
2  94 

295 

296 

297 
293 

O  0  <5 


312 
3  1  3 
3  1  4 
3  15 
9  i  6 
1  7 


C 

C 


c 

c 

c 

c 


than  or  equal  to  the  square  root  of  Const. 

This  is  found  by  finding  the  Magnitude  of  the  hypotneuse  of  a 
right  triangle  with  base  and  height  of  the  current  row 
and  colunn .  The  slope  is  the  row  (  Ikount  )  /  coluMn  (  K  ). 
This  procedure  is  followed  until  the  upper  bound  is  reached. 


C  *  *  * *  *  *  ****** 

C 

100  Do  140  Ikount  =  16 
Do  130  K  =  0 
If  <  K  .gt 


-  1 


16 


Ikount  )  Goto  140 
Run  =  FLOAT  (  K  ) 

Rise  =  FLOAT  <  Ikount  ) 


C* *************************************  ********************************* 

c 

C  Get  angle  for  this  slope  and  be  sure  it  does  not  fall  outside 

C  the  bounds  on  theta. 

C 

C**t***t* ****************************************  *********************** 

c 

Angle  =  ATAN  (  Run  /  Rise  ) 

If  <  Angle  It.  Radlow  )  Goto  130 
If  <  Angle  .gt.  Radhi  >  Goto  140 


C* *************  *****************************  **************************** 

c 

Deter  Mine  how  Many  coluMns  out  to  "run"  (  Jinc  )  and  how  Many 
rows  to  "rise"  (  Ikount  ).  Center  the  circle  on  the  point 
123,123  of  RSourc  and  place  a  one  at  locations  where  pupil  is 
to  transmit. 


C* **********  ************ 

C 


*  *  *  *  ************************ 


30  0 

Jinc  =  K 

30  1 

Islope  =  Ikount 

3  02 

J  r  o  w  =  0 

303 

Do  120  I  =  0  ,  16  ,  Islope 

304 

Rlong  =  I  *  I  +  Jrow  *  Jr o  w 

305 

If  (  Rlong  gt.  Const  )  Goto 

130 

306 

If  (  Rlong  It.  Other  )  Goto  115 

307  ' 

RApert  (  128  -  I  ,  128  +  Jrow  )  = 

1  .  0 

3  0  3 

1  15 

Jrow  =  Jrow  +  Jinc 

309 

120 

Continue 

3  1  0 

130 

Continue 

3  1  1 

140 

Continue 

c**t*  *  ****************************************************************** 

c 

C  Do  the  sane  as  above  for  the  rest  of  the  pupil. 


u 


C  *  *  t  *  *  *  *  **********************************  ******  ************************ 
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319 

150 

Do  180  Jkount  =  16 

.  1  ,  - 

i 

320 

Do  170  K  =  16  ,  - 

16  .  - 

1 

.321 

If  <  ABS  (  K  ) 

.gt.  Jkount  )  Goto 

322 

Runovr  =  FLOAT 

<  K  ) 

323 

Riseov  =  FLOAT 

(  Jkount 

) 

324 

If  (  K  .  ge .  0  ) 

Then 

325 

Angle  =  (  Pi  / 

2 .  >  - 

ATAN  (  Run 

326 

Else 

327 

Angle  *  (  Pi  / 

2.  )  - 

ATAN  <  Run 

323 

Endif 

329 

If  <  Angle  .It. 

R  a  d  1  o  w 

)  Goto  170 

330 

If  <  Angle  .  g  t  . 

Radhi  ) 

Goto  180 

331 

If  (  K  .eq .  0  ) 

Then 

332 

Slope  *  0  . 

33  3 

Else 

3  3  4 

Slope  =  Riseov  /  Runovr 

335 

Endif 

336 

If  (  Slope  .  e q  . 

0  )  Then 

337 

Do  155  J  =  0  , 

16 

333 

J2  =  J  *  J 

339 

If  (  <  J2  It 

Other 

)  or  (  J 

340 

RApert  (  123 

,  128  + 

J  >  =  i  .  0 

341 

155 

Continue 

342 

Else 

343 

Ir  ow  =  0 

344 

line  =  -  K 

345 

Do  160  J  =  0  , 

16  ,  Jkount 

346 

Rlong  =  J  *  J 

+  Irow  * 

Irow 

347 

If  (  Rlong  gt 

Const 

)  Goto  170 

348 

If  (  Rlong  It 

Other 

)  Goto  156 

3  49 

RApert  <  128  + 

Irow  , 

128  +  J  >  : 

350 

156 

Irow  =  Irow  + 

line 

351 

16  0 

Continue 

352 

Endif 

353 

170 

Continue 

35  4 

130 

Continue 

355 

190 

Do  230  Ikount  «  16 

,  1  ,  - 

1 

356 

Do  220  K  =  16  ,  0 

,  -  i 

357 

If  <  Ikount  .It 

.  K  )  Goto  220 

358 

Rise  =  FLOAT  ( 

Ikount  ) 

359 

Pun  =  FLOAT  (  K 

) 

360 

Angle  =  Pi  -  ATAN  (  Run 

/  Rise  ) 

361 

If  (  Angle  .it 

P  a  d  1  o  w 

)  Goto  220 

362 

If  (  Angle  gt 

Radhi  ) 

Goto  230 

363 

If  <  Run  eq.  0 

)  Pun  = 

0  .  1 

364 

Slope  =  -  Rise 

/  Run 

365 

If  (  Slope  It. 

-  16  ) 

Then 

3  6  6 

Islope  =  -  1 

3  6  7 

Jinc  =  0 

368 

Jrou  =  0 

369 

Else 

3  T1  0 

Islope  =  Ikount 

7  1 

Jinc  =  K 

170 


/  Riseov  ) 


/  Riseov  ) 


g  t  .  Const  )  )  Goto  155 


1  .  0 
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72  Jrow  =  0 

73  Endif 

74  Do  210  1=0  ,  -  16  ,  -  Islope 

75  Rlong=I*I+  Jrow*  Jrow 

76  If  (  Rlong  .gt.  Const  )  Goto  220 

77  If  <  Rlong  .It.  Other  )  Goto  205 

73  RApert  <  128  -  I  ,  128  +  Jrow  )  =  1.0 

7?  205  Jrow  =  Jrow  +  Jinc 

30  210  Continue 

31  220  Continue 

82  230  Continue 

93  C 

3  4  CMM  ******************  ************************************************* 

35  C 

86  C  Since  pupil  is  s  y  n Metric,  fill  in  other  half  just  like  the  other 

87  C 

33  CM* ***********************************************  ********************* 

39  C  1 

70  240  Do  260  J=  16  ,  1  ,  -1 

91  Do  250  1=0  ,16 

92  RApert  <  128  +  I  ,  128  -  J  >  =  RApert  (  128  +  I  ,  128  +  J  ) 

93  RApert  <  128  -  I  ,  123  -  J  >  =  RApert  (  128  -  I  ,  128  +  J  ) 

94  250  Continue 

95  260  Continue 

96  C 

97  C***«*» *************************** t ********************  ***************** 

9  3  C 

99  C  See  if  operator  would  like  to  add  to  the  aperture. 

0  0  C 

01  C  **********************************************  ************************* 

02  C 

03  Write  <  7  ,  261  ) 

04  261  For nat(/,$, 'Would  you  like  to  add  to  the  aperture  [  y/n  ]  ?  '  ) 

05  Read  (  5  ,  262  )  Answer 

06  2  &  2  Forriat(Al) 

07  If  (  Answer  .eq.  y  )  Goto  1 

0  3  C 

09  C  **********************  ***************************  ********************** 

1 0  c 

11  C  See  if  a  picture  is  wanted. 

1  2  C 

13  CMMMMMMMMMMMMMMMMM  **********  MM*  MM*  ****************  * 

1  4  C 


Write  <  7  ,  270  ) 


1  to 

2  7  0 

Fornat(/,*, 

'Do 

you  want 

t  h 

e  aperture  func 

tion  plotted  on  the 

scree 

1  7 

In  t  y/n  ]  7 

'  ) 

13 

Read  (  5  , 

230 

)  Answer 

19 

23  0 

Foriiat  (A1  ) 

20 

If  (  Answer 

.  n  e 

y  )  Goto 

290 

2  1 

Call  P  1  tr  s  t 

(  2 

,  5  ) 

7} 

il  — 

290 

Write  <  7  , 

30  0 

) 

2  3 

30  0 

Forna  t ( /  ,  J , 

'Do 

you  want 

a 

hardcopy  of  the 

aperture  function  [ 

y/n 
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1 

2 

3 

4 

5 

6 
/ 
8 
? 

10 
1 1 
12 

13 

14 

15 

16 
17 
13 
1? 
20 
21 
t—  w 

23 

24 

25 

26 
27 
23 

2  9 

3  0 

31 

32 

33 

34 

3  5 

36 

37 
2  9 
39 

4  0 

41 

42 


Subroutine  Invert  (  Ichose  ) 


C 


c 

C  This  subroutine  per -for  ms  Magic  on  the  source  array 

C  enabling  the  low  frequency  components  to  be  displayed  in  the 

C  center  as  opposed  to  the  ends  as  would  otherwise  be  the  case. 

C 

C** ********************************  ************************************* 


c 

Common  /  A  r  g  s  /  Real  (  256  )  ,  Rimag  <  256  )  ,  RSourc  (  256  ,  25 

16  )  ,  CSourc  <  256  ,  256  )  ,  RApert  <  256  ,  256  )  , 

2  Plot  (  256  ,  256  )  ,  Ar  (  2  )  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 


C********** 


C 

C  Every  other  element  of  the  arrays  RSourc  and  CSourc 

C  starting  with  the  second  element  is  multiplied  by  -1. 

C  The  preprogrammed  sources  (  point,  two  point,  edge,  slit  )  do 

C  not  need  to  be  inverted  since  their  transforms  were  calculated 

C  analytically  which  put  the  low  frequencies  in  the  middle. 


C 

C  ****************************************  ******************************* 

c 

If  (  Ichose  .It.  5  )  Goto  50 
W  r  i  t  e  (  7 , 3  ) 

3  For mat</, 'Inverting.') 

Do  20  J  =  1  ,  256  ,  2 

Do  10  I  =  2  ,  256  ,  2 

RSourc  <  I  ,  J  )  =  -RSourc  (  I  ,  J  ) 

CSourc  <  I  ,  J  )  -  -  CSourc  <  I  ,  J  ) 

10  Continue 

20  Continue 
C 


Do  40  J  =  2  ,  256  ,2 

Do  30  I  =  1  ,  256  ,  2 

RSourc  <  I  ,  J  ) 
CSourc  <  I  ,  J  ) 


RSourc  (  I  ,  J  ) 
CSourc  <  I  ,  J  ) 


30  Continue 

40  Continue 
50  Return 
End 


0  Errors  detected 
42  Source  lines  read 
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u  M  m  M  u  *  m  »  M  U  M  M  MJ  M  »  M  » 


This  prograw  takes  the  2-Dimensional  FFT  of  the  source 
ir radiance  distribution.  This  is  done  analytically  for  point, 
two  point,  edge,  and  slit  sources.  Other  sources  are  done 
using  a  one-dimensional  fast  Fourier  transform  (FFT)  program. 
The  rows  of  the  source  are  transformed  first  3nd  the  columns 
second.  The  transform  may  be  displayed  on  the  terminal  and 
plotted  if  desired. 


C  *  *  * 

c 


Common  /  Args  /  Real  <  256  )  ,  Rimag  (  256  )  ,  RSourc  (  256  ,  25 

16  )  ,  CSourc  <  256  ,  256  )  ,  RApert  (  256  ,  256  )  , 

2  Plot  (  256  ,  256  )  ,  Ar  (  2  )  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 
C 

Character  Answer  ,  y 
Data  y/'y'/ 

C 

C  *  m  H  t  M  M  m  M  t  M  M  U  U  M  t  M  i  I  M  M  M  H  1 1  M  1 1  i  U  M  1 1  U  1 1  t  m  t  M  1 1  tl  >  M 

C 

C  Let  the  user  know  that  fftsrc  has  been  invoked. 

C 

C  M  rn  I  t  ti  M  i  M  U  U  M  M  k  *  M  M  M  M  M  M  *  M  I  H  t  t  M  •  1 1  ti  M  I  i  M  M  I  I  1 1 1  t  (  t  t  i  I  I 

C 

Write  (7,5) 

5  For  n  a  t  < /  ,  'Transforming  source.') 


C 

c  *  *  * 

c 


C  Ichose  is  passed  into  fftsrc  from  the  calling  program  (synapt). 

C  The  value  of  Ichose  corresponds  to  a  particular  type  of  source. 

C 

C  Ichose  =  1  Point  Source 

1  Ichose  =  2  Two  Point  Source 

L  Ichose  =  3  Edge  Source 

C  Ichose  =  4  Slit  Source 

L  Ichose  =  5  Circular  Source 

C  Ichose  =  6  Anything  Else 

C 

C  The  computed  Goto  statement  transfers  the  program  to  the 

C  appropriate  part  of  the  program  to  calculate  that  source's 

C  transform. 

r- 

u 

C*tt ********************************************  ************************ 

L 

Goto  (  10  ,  40  ,  30  ,  80  ,  490  ,  490  )  Ichose 

C 

c  **********  **********************************************  *************** 

c 

C  Point  source  has  a  plane  wave  as  its  transform.  Make  RSourc 
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equal  to  1  at  all  locations  and  CSourc  equal  to  0  (  phase 


Do  30  J  =  1  ,  256 


59 

Do  20  I 

S 

1 

,  256 

60 

RSourc 

( 

I 

,  J  )  =  1  .  0 

6  1 

CSourc 

( 

I 

,  J  >  =  0  .  0 

62 

2  0 

Continue 

63 

30 

Continue 

64 

Goto  560 

C 

C  *  «  *  : 

C 

c 

c 

c 

c 

c 

c  *  *  *  ■ 
c 

40 


Two  point  source  transform:  is  a  cosine  wave.  The  two  points 
are  separated  by  a  distance  specified  by  the  user.  Therefore, 
the  cosine  wave  is  a  Cos  <  N  *  pi  /  (  256  /  Separation  )  ). 

The  columns  of  RSourc  are  set  equal  to  twice  this  quantity. 

Xstart  =  -  (  Float  <  Length  )  /  2.  )  *  3.1415926535 

Xconst  =  3.1415926535  /  (  256.  /  Float  <  Length  )  ) 

Do  50  I  =  1  ,  256 
X  =  Xstart  +  Float  <  I  )  *  Xconst 


Real  <  I  )  = 
Continue 
Do  70  J  =  1 
Do  60  I  =  1 
RSourc  (  I 
Continue 
Continue 
Goto  560 


t  Cos  <  X  ) 


256 
256 
,  J  ) 


Real  (  I  ) 


C  *  M 

c 


******************************** 


Fourier  transforms  of  an  edge  or  a  sli 
has  a  finite  width  (  I  wide  )  and  an  in 
This  slit  has  both  finite  width  and  le 
edge  is  a  Sine  (  pi  *  X  /  (  256  /  I w i d 
periodicity  of  the  edge.  The  transfer 
calculated  except  that  it  is  the  produ 
The  edge  is  oriented  such  that  the  Ion 
azimuth  and  its  width  along  the  range, 
the  same  way. 


of  an  edge  or  a  slit  are  similar.  The  edge 
(  I  wide  )  and  an  infinite  length  (  Length  ). 
finite  width  and  length.  The  transform  of  an 
*  X  /  (  256  /  I  wide  )  )  where  256  is  the 
edge.  The  transform  of  the  slit  is  similarly 
hat  it  is  the  product  of  two  sine  functions, 
d  such  that  the  long  side  runs  along  the 
th  alona  the  range.  The  slit  is  oriented 


C  *  *  ♦ 

C 

SO 


Y const  =  Float  (  256  /  Iwide  ) 
xconst  =  Float  (  256  /  Length 
rstart  =  -  128.  *  3.1415926535 

Xstart  =  -  128.  *  3  1415926535 
Xincr  =  3.1415926535  /  Xconst 


/  Yconst 
/  Xconst 
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107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 
12? 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 
14  7 
148 
14? 

150 

151 

152 

153 

154 

155 

156 

157 
153 
15? 


Yincr  =  3.1415926535  /  Yconst 
Do  90  I  =  1  ,  256 

Y i  =  Ystart  +  Float  (  I  -  1  )  *  Yincr 
If  (  Abs  (  Yi  )  gt.  IE-7  )  Then 
Sine  =  (  Sin  (  Yi  )  /  Yi  ) 

Else 

Sine  =  1.0 
En d  i  f 

Real  (  I  )  =  Float  <  Length  >  »  Sine 
90  Continue 

Do  110  J  =  1  ,  256 

Do  100  1=1,  256 
RSourc  (  I  ,  J  )  =  Real  (  J  ) 
lOOContinue 
110  Continue 

Do  120  1=1,  256 

X  =  Xstart  +  Float  (  I  -  1  )  »  Xincr 
If  (  Abs  (  X  )  gt.  IE-7  )  Then 
Sinc=(Sin(X)/X) 

Else 

Sine  =  1.0 
End  if 

Real  (  I  )  =  Float  (  Iwide  )  *  Sine 
120  Continue 

Do  140  J  =  1  ,  256 
Do  130  I  =  1  ,  256 

RSourc  (  I  ,  J  )  =  RSourc  (  I  ,  J  )  *  Real  (  I  ) 

130  Continue 
140  Continue 
Goto  560 
C 

c  |  m  nn  t  t  M  M  M  »  M  t  t  t  M  t  M  M  »  M  M  ♦  M  *  M  *  M  I  H  U  M  t  t  M  t  t  M  <  I  I  M  I  I  I  U  M  I 

c 

C  For  anything  else,  calculate  the  FFT  First  the  rows. 

C 

Cm  Ut  H  t  ♦  tM  U  M  M  »  M  M  *  M  M  MM  t  HUM  M  H  t  J  MM  M  m  M  m  M  I  M  M  M  M  I 

C 

490  Do  520  1=1,  256 

Do  500  J  =  1  ,  256 

Real  (  J  >  =  RSourc  (  I  ,  J  ) 

Rinag  (  J  )  =  CSourc  (  I  ,  J  ) 

500  Continue 

Call  FFT  (  Real  ,  Rinag  ,  256  ,  i  ) 

Do  510  J  =  1  ,  256 

RSourc  (  I  ,  J  )  =  Real  <  J  ) 

CSourc  (  I  ,  J  )  =  Rifiag  (  J  ) 

510  Continue 

520  Continue 
C 

CfM*Mt*t*tM*ttt***tt*tt*t*t«ttt**t*ttt*****«*ttttt*t*t**M*ttttttt*tt 

C 

C  Now  the  colunns. 

C 
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60 

C  *  *  * 

61 

C 

62 

63 

64 

t>5 

t>6 

530 

o7 

68 

o9 

70 

71 

540 

72 

550 

73 

c 

74 

ClM 

75 

c 

76 

c 

77 

c 

78 

C  t  *  1 

79 

c 

80 

560 

31 

570 

82 

83 

84 

580 

85 

86 

87 

5  9  0 

98 

600 

89 

9  0 

91 

92 

93 

94 

95 

Do  550  J  =  1  ,  256 

Do  530  1=1,  256 

Real  (  I  )  =  R  S  o  u  r  c  <  I  ,  J  > 
Rifiag  <  I  )  =  CSourc  (  I  ,  J  ) 

Continue 

Call  FFT  (  Real  ,  RiMag  ,  256  ,  1  ) 

Do  540  1=1,  256 

RSourc  (  I  ,  J  )  =  Real  (  I  ) 
CSourc  (  I  ,  J  )  =  RiMag  (  I  > 

Continue 

Continue 


See  if  the  transform  is  to  be  displayed  or  plotted. 

t  M  t  M  *  M  M  H  t  H  M  U  I  M  M  M  1 1  It  U  t  *  I  1 1 1  1 1 1 1 1 1 1 1 1 


Write  (  7  ,  570  ) 

For  M3 1 (  /  ,1  ,  'Do  you  want  the  source  FFT  displayed  on  the  screen  [  y 
1  /  n  1  ?  '  ) 

P  e  a  d  (  5  ,  580  )  Answer 
FOf'M3t<Al  ) 

If  (  Answer  .  n  e  .  y  )  Goto  590 
Call  P  1 1 r s t (  3  ,  0  ) 

Write  (  7  ,  600  ) 

Format  (/,$,'  D  o  you  want  a  hardcopy  of  the  source  FFT  t  y/n  1  ?  ') 

Read  (  5  ,  580  )  Answer 

If  (  Answer  .ne.  y  )  Return 

Call  Plotin 

Call  Fltrst  (  3  ,  0  ) 

Call  Flotof 

Return 

End 


0  Errors  detected 
195  Source  lines  read 
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i 


i 


u 


> 


8 

9 

10 
1  i 
12 
13 

1  4 
15 
1  6 

17 

18 
1? 

2  0 
21 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3  4 
35 
35 

t  7 

3  8 
3-5 

40 

41 

4  2 


C 

c  *  *  * 

c 

c 

c 

c 

c 

c 

c 

C  *  *  * 

c 


Subroutine  IFTSrc 
**************************** 


This  subroutine  calculates  the  inverse  Fourier  transform  of  the 
product  of  the  aperture  and  source  transforms.  This  product  is 
stored  in  RSourc  and  CSourc.  This  inverse  transform  is  the 
goal  of  this  program.  The  inverse  transform  can  be  displayed  on 
the  screen  or  a  hardcopy  output  may  be  obtained. 


******* 


****************************** 

25 


Common  /  Ar  g  s  /  Real  (  256  )  ,  Rimag  .(  256  )  ,  RSourc  (  256 

16  )  ,  CSourc  (  256  ,  256  )  ,  PApert  (  256  ,  256  )  , 

2  Plot  (  256  ,  256  )  ,  Ar  (  2  )  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 


Character  Answer 
Data  y / ' y ' / 


29 


10 


2  0 
30 


Hr  i  te (7 , 29 ) 

Formate/, 'Inverse  transforming.  '  > 

Do  30  I  =  1  ,  256 

Do  10  J  =  1  ,  256 

Real  <  J  )  =  RSourc  (  I  ,  J  ) 
Rimag  (  J  )  =  CSourc  (  I  ,  J  > 

Continue 

Call  FFT  (  Real  ,  Rimag  ,  256  ,  -  1  ) 

Do  20  J  =  1  ,  256 

RSourc  (  I  ,  J  )  =  Real  (  J  ) 
CSourc  <  I  ,  J  )  =  Rimag  (  J  ) 

Conti.  •>  e 
Continue 


J  ) 

J  > 


4  0 


Do  60  J  =  1  ,  256 

Do  40  I  =  1  ,  256 

Real  (  I  )  =  RSourc  <  I 
Rimag  (  I  )  =  CSourc  (  I 

Continue 

Call  FFT  (  Real  ,  Rimag  ,  256  ,  -  1  ) 

Do  50  I  =  1  ,  256 

RSourc  <  I  ,  J  )  =  Rea  1  (  I  > 
CSourc  <  I  ,  J  )  =  Rimag  (  I  ) 


43 

50 

Continue 

44 

6  0 

Continue 

45 

1 1  r  i  t  e  <  7  , 

70  > 

46 

70 

Fornat(/,l, 

'Do  you  want  the  source 

47 

1  /  n  ]  ?  '  ) 

• 

48 

Read  (  5  , 

80  )  Answer 

,  - 

4? 

80 

Fornat(Al) 

L  , 

W* 

50 

If  (  Answer 

.ne.  y  )  Goto  90 

w*. 

51 

Call  Pltrst 

(5,0) 

V  ’ 

52 

9  0 

Write  (  7  , 

100  ) 

•  .  53 

i  e  o 

Format(/, t , 

'Do  you  want  a  hardcopy 

• 

r 

if 


Fortran/9000 
Ve  r .  4.02 


Mon  Sep  23  09:45:31  1985 
f  f  t  .  f  Page  1 


Subroutine  FFT  (  FP. 


IDIR  ) 


Data  is  in  FR  (real)  and  FI  (imaginary)  arrays. 
Computation  is  in  place,  output  replaces  input. 
Number  of  points  must  be  N. 

FR  and  FI  must  be  dimensioned  in  the  main  program 
IDIR  =  +  1  (  spatial  to  frequency  ) 

IDIR  =  -  1  <  frequency  to  spatial  ) 


Dimension  FR ( 1 )  ,  FI(1) 

If  (  IDIR  .  EQ  .  +  1  )  goto  10 

Do  9  JJ  =  1  ,  N 

FI  <  JJ  )  =  -1 .  *  FI  ( JJ  ) 


FI  ( JJ  )  =  -1  .  *  FI  ( JJ  ) 

MR  =  0 

UN  =  N  -  1 

Do  2  M  =  1  ,  NN 

L  =  N 

L  =  L  /  2 

If  (  MR  +  L  .  GT.  NN  )  goto  1 
MR  =  Hod  (  MR  ,  L  )  +  L 
If  <  HR  . LE.  M)  goto  2 
TR  =  FR  <  M  +  1  > 

FR  (  M  +  1  )  =  FR  (  MR  +  1  ) 

FR  (  MR  +  1  )  =  TR 
TI  =  FI  <  M  +  1  > 

FI  <  M  +  1  )  =  FI  (  MR  +  1  ) 

FI  <  NR  +  1  )  =  TI 
Continue 
L  =  1 

If  (  L  .GE.  N  )  goto  12 
ISTEP  =  2  *  L 


32 

EL 

= 

L 

3  3 

Do 

4 

M  =  1  ,  L 

34 

A 

=  3 

. 1415926535 

♦Float  <  1  -  h  >  / 

35 

MR 

= 

Cos  <  A  ) 

36 

HI 

= 

Sin  (  A  > 

37 

D  o 

4 

I  =  M  ,  N 

f 

ISTEP 

3  8 

J 

=  I 
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1  Subroutine  Grph  of 

2  C 


3 

4 

5 
A 


C* **********************  **************** 


c 

C  This  subroutine  turns  off  the  graphics  when  the  Main  progran 

C  is  finished. 


8  C  *  *  *************************  ******************************************** 

9  C 

10  Co  mm  on  /  Args  /  Real  (  256  )  ,  RiMag  (  256  )  ,  R.Sourc  <  256  ,  25 

11  lb  )  ,  CSourc  (  256  ,  256  )  ,  RApert  (  256  ,  256  )  , 

12  2  Plot  <  256  ,  256  )  ,  Ar  (  2  )  ,  Lower  ,  Iupper  ,  Length  ,  Iwide 

13  C 

14  Call  JWOFF  <  1  ) 

15  Call  JWEND  <  1  ) 

16  Call  JEND 

17  Return 

13  End 


0  Errors  detected 
18  Source  lines  read 


The  following  pages  are  taken  directly  from  the  HP  2623a  graphics 
terminal  through  the  use  of  its  internal  printer.  Two  sample  runs  are 
illustrated.  The  first  is  a  flawless  execution  of  the  program  to 
illustrate  what  happens  when  everything  goes  right.  The  second  sample 
shows  what  can  happen  when  everything  goes  wrong  illustrating  the  vari¬ 
ous  error  messages  and  what  the  model  does  to  recover. 

Anyone  wishing  to  run  the  model  will  first  have  to  visit  the  system 
manager  of  the  HP  9000  computer.  He  is  Jeff  Sweet  and  is  located  in 
Bldg  622  along  with  the  computer  and  has  an  office  in  room  104.  He  can 
be  reached  by  phone  at  extension  5-6361.  He  will  assist  the  user  in 
getting  an  account  on  the  computer  and  getting  the  user  familiarized 
with  the  system  in  general.  Once  logged  in  and  in  the  proper  account  as 
per  his  instructions,  the  model  can  be  invoked  by  typing  'synapt' 
followed  by  a  carriage  return.  The  sample  model  runs  begin  at  this 
point . 


i  /users/ kane/graph  i .  ■  ^  u  j _ p/iiap  i  .:?•/**•  .-1*1  i  ‘  -  •  i*. 

Enter  your  source  irradiance  distribution.  You  may  choose  from  one  of i'.the' pre¬ 
programmed  distributions  below  or  create  your  own.  Type  the  number  of  j  your;,*j '* ’ *;* 
selection  after  finding  your  choice  on  the  menu  below  when  prompted ..  >•  V.’T  / 

point  source:  I  .  •  •  i 

h  two  point  source:  2  "i  '■  '  :  v  -  i  ■  •  :  u*’ 

fln  edge:  3  :  ’■  '  ■  -j-  -  ;V  - 

s  siit-  .  ;  d  -m:  :■ 

h  circle  of  variable  radius:  5  ’  '  '  V ' i  f  • 

Your  own-  creation:  6  I  .  ■  •  •  ••  f7;  •  'r- 1 :1\ •!..*?  •’>. 

n  ■'  ■  ■  ■,  ;  ••  '•  *  • ;  |.;-i  ,.}■*!  i 

Enter  your  selection  [1-6]  :  1  J  . * ? ii'.-’ v|;  : ’ 


mapt 


point  source:  1  ' 

h  two  point  source:  2 

An  edge :  3 

H  slit:  4 

fl  circle  of  variable  radius:  5 

Your  own • creation :  6 

Enter  your  selection  [1-6]:  1 

Du  you  want  the  source  plotted  on  the  screen  1  y/n  ]  ?  n 
Do  you  want  a  hardcopy  of  the  results  [  y/n  1  ?  n 


enter  range  (Km)  and  lens  separation  <m)  as  x.x  x.xx  :  1.0  0.50 
Enter  collector  speed  (ft/s)  and  stability  as  xxxx.x  xx.x  :  Q88C 
Enter  theta  for  your  particular  aperture. 

What  range  does  theta  lie  within  from  0  to  180  degrees?  • 

[inter  your  answer  as  nnn  nnn  :  037  143  ' — ’-’T' 

would  you  like  to  add  to  the  aperture  [  y/n  ]  ?  n  i 


■ViWMfrr.J-i 

■  I  .  ;  • ’  .  I  i  '  ,-i  *  ! 

'  h  •;!.>  |  ■ 

:a/! 

AS,  .  f.v, 


lu  you  want  the  aperture  function  plotted  on  the  screen  [  y/n  ],.  ?  n 

'  *  f  :■  .  i{ 

Jf  \  '  >  ;  *«. 

Du  you  want  a  hardcopy  of  the  aperture  function  L  y/n  ]  ?  n  - 


Iranstorming  source.  ■■■ 

Do  you  want  the  source  EFT  displayed  on  the  screen  [  y/n  ]  ?:,n 
Do  you  want  a  hardcopy  of  the  source  FFT  [  y/n  ]  ?  n 
Multiplying  pupil  by  source  FFT.  ; 


0880.0  10.0 -j&v  :A 
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;  l.; M ,.:*H*-* "V;’ 
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y  ■ .  -  •  •  ’.;i.  .*jv 

v’  !i  '  ’  :  * 


■  •  •  1  \  r  '■  4fV  v 


Do  you  want  a  plot  of  the  product  of  the  FFT  of  the  source  and  ,  the  aperture  .f. . 
iJistributions  [  y/n  ]  7  n  :  ■  1 

Do  you  want  a  hardcopy  of  the  product  [  y/n  ]  ?  n  \  ,  ■  E1  ; 1 

\  .. 

Inverting.  ^  '"’;7 

>  .  V  ■r^X^r. 

Inverse  transforming.  .  •  1  f  v  l,* 

■"  l\;  •  r" 

Jo  you  want  the  source  IFT  displayed  on  the  screen  t  y/n  J  ?  n  f**;. 

Do  you  want  a  hardcopy  of  the  source  IFT  [  y/n  J  .?  n  :  '  i  -  !■->■■■ 

>  .  i  ,  l;  •*,,  */ 

,  ;  '  <  .  ■  .  '  *  jj  4,  f,t  J  j 

i*  «you  desire  to  try  another  source  and  aperture  [  y/n  ].  ?  n'  if <iJ.<  D 
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users/kane/graph 


•v>r1 


7user s/kane/granh  tRO'i  -  Isynapt 


a 


Enter  your  source  irradiance  distribution.  You  may  choose  from  one  of  the  pre- 
.yogrammed  distributions  below  or  create  your  own.  Type  the  number  of  your 
s-ilection  after  finding  your  choice  on  the  menu  below  when  prompted. 


P  point  source: 

P  two  point  source: 

.On  edge: 

0  slit: 

n  circle  of  variablle  radius 
Your  own  creation: 


1 

O 

L. 

3 

4 


enter  your  selection  11-61:  1 


Oo  you  want  the  source  plotted  on  the  screen  L  y/n  1  ?  n 


bo  you  want  a  hardcopy  of  the  results  1  y/n  j  ?  n 

Enter  range  (Km)  and  lens  separation  <m)  as  x.x  x.xx  :  1.,.5 

Enter  collector  speed  (ft/s)  and  stability  as  xxxx.x  xx.x  :  880. ,10. 


► 

r. 

k 


Enter  theta  for  your  particular  aperture. 

what  range  does  theta  lie  within  from  0  to  180  degrees? 

f..w.ter  your  answer  as  nnn  nun  :  037  143 

You  have  exceeded  the  maximum  Dtheta  of  .008  degrees  for  the  range  entered. 
Try  again. 


Enter  range  (Km)  and  lens  separation  (in)  as  x.x  x.xx  :  1.00  1  .000 

jtav  Enter'  collector  speed  (ft/s)  and  stability  as  xxxx.x  xx.x  :  888888  1 

Enter  theta  tor  your  particular  aperture.  *. 

Uhat  range  does  theta  lie  within  from  U  to  180  degrees? 

“  Enter  your  answer  as  nnn  nnn  :  37  180 
^  Your  combination  of  lens  separation  and  range  caused 

the  upper  spatial  frequency  to  fall  outside  the  aperture. 

E.  Please  try  again. 

Enter  range  (Km)  and  lens  separation  (m)  as  x.x  x.xx  : 
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// 
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This  thesis  was  concerned  with  the  development  of  a  computer  model 
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The  system  modeled  consisted  of  two  lenses  mounted  on  a  movable  platform. 

The  lenses  were  separated  by  a  fixed  distance  and  travelled  in  a  direction 
parallel  to  this  separation.  The  coherence  of  radiation  present  at  each  lens  emana¬ 
ting  from  a  real  source  was  measured  yielding  the  Fourier  transform  of  the  source 
intensity  distribution,  according  to  the  van  Cittert-Zernike  theorem  (2:510). 

The  transform  was  then  multiplied  by  an  effective  aperture  (obtained  from  the  motion 
and  position  of  the  lenses  relative  to  the  source).  An  inverse  Fourier  transform  was 
then  applied  to  this  result  yielding  the  image.  This  is  the  process  modeled  by 
the  computer. 

The  results  indicated  that  new  means  of  image  interpretation  must  be  developed 
in  order  to  make  the  results  useful.  This  is  due  to  the  fact  that  the  system  behaves 
much  like  a  high  pass  filter  and  the  image  is  edge  enhanced  and  not  a  scaled 
version  of  the  geometric  image. 
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